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Remote Sensing offers promise for studying sea ice conditions and the processes 
which change them. At present, there is too little common ground between remotely 
observable data and the quantities described in sea ice models. Our interest is in 
identifyimng measurements which can be made from satellites and which describe 
properties relevant to the balances of mass, heat and momentum of the ice cover. As a 
prerequisite, we wish to describe the random spatial structure of sea ice and its field of 
motion, so that adequate sampling strategies can be devised. 

The research performed under this grant has been focussed on the piece-like 
structure of the ice pack, as described by the distribution of floe sizes. Objectives were 1} 
to clarify how the several useful definitions of floe size distribution are interrelated, 2) to 
consider the practicality of different measurement techniques, 3) to make measurements 
of typical distributions, and 4) to investigate the effect of sample size on the sampling 
error. 

These objectives have been accomplished. The results are contained in a paper 
daringly entitled "Sea Ice Floe Size Distribution" and attached here as Appendix I. 
Summertime distributions of N(p), the number of floes per unit area with diameter no 
smaller than p, behaves roughly as p“ with a ranging from -1.7 to -2.5. The sampling 
theory gives the variance of an estimate for the fraction g of a test area with diameter D 
covered by floes of diameter d as g(l—g)N~ l , where Nv^{D/d) z can be thought of as the 
number of independent samples in the test area. This relationship and similar ones for 
other test regions have been copnfirmed by measurements of sea ice images: LANDSAT 
visual images, U-2 aerial photgraphs, and a mosaic of NASA aerial photographs. 

The basic procedures and theoretical facts are in hand for measuring floe size 
distribution. Different distributions can be defined and interrelated. Distributions have 
been measured for several stages of break-up-all of which show distinct floes. The 
extension to winter conditions where floes are not well defined should be pursued; a useful 
approach would be a study of lead geometries. Another useful step would be the 
application of these measurement techniques to many images of sea ice to observe the 
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seasonal and regional variability of floe size distribution. 

The sampling theory devised for floe size distribution has application to any problem 
involving spatial sampling. It was applied to the sampling of ice thickness distribution, as 
expounded in the paper "Ice Thickness Distribution— Theory and Measurement" attached 
as Appendix H. The principle result is that for thick ice, samples further apart than about 
a hundred meters are independent, while independent samples of thin ice (thinner than 
one meter) must be over a kilometer apart. 

The velocity field of sea ice is demonstrably the movement of rigid pieces. The 
consequences of this on the spatial statistics of the field are summarized in the paper 
"Kinematics of Sea Ice" attached as Appendix IH. Included is the description of a model of 
the piecewise rigid body motion and its effect on the estimation of open water formation 
and ridging from a “large scale" measurement of deformation. The pieces are taken to be 
defined by a random set of lines (a Poisson field) with an average piece size of 16 km, 
taken from synthetic aperture radar data of Hall and Rothrock (1981). The differential 
movement of the pieces is given a Gaussian distribution. Many realizations of these pieces 
and their motion are simulated. For each, the total opening and closing at the floe 
boundaries is plotted versus a "large scale" estimate of deformation made from the 
velocity at three points. Each realization gives one data point. The scatter in the data 
points is several hundred percent and shows how poorly a three point strain rate estimate 
allows us to determine the actual small scale opening and closing. 

RECOMMENDATIONS FOR FUTURE WORK 

Floe size measurements should be extended into winter. The most attractive 
approaches are the study of the branching geometries of leads and the identification of 
pieces through their motion rather than their appearance. It would be useful, at this 
stage, to measure floe size distributions from a variety of ice conditions, and to begin to 
observe the spatial and temporal variability of this property of sea ice. 
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Further work, currently in progress, is aimed at direct observation of local (small 
scale) opening and closing from high resolution velocity field obtained from SAR imagery. 
SAR can also provide a better understanding of the statistics of the rigid body motion of 
pieces: for instance, whether their rotation and translation depends on floe size. 

REFERENCE 

Hall, R.T. and D.A Rothrock (1981) Sea Ice Displacements from Synthetic Aperture Radar, 
Journal of Geophys. Res,, 86, 11078-82. 

PUBLICATIONS ISSUED DURING THIS GRANT 

D. A, Rothrock and A. S. Thorndike, Sea Ice Floe Size Distribution, submitted to 
Journal of Geophys. Res., attached as Appendix I. 

D. A Rothrock, Ice Thickness Distribution— Measurement and Theory. NATO Advanced 
Study Institute, Maratea, Italy; September 1981. In The Geophysics of Sea Ice, N. 
Untersteiner, editor. (In preparation, Plenum Publishing, New York.) Attached as 
Appendix II, 

A S. Thorndike, Kinematics of Sea Ice. NATO Advanced Study Institute, Maratea, 
Italy; September 19B1. In The Geophysics of Sea Ice, N. Untersteiner, editor. (In 
preparation. Plenum Publishing, New York.) Attached as Appendix El. 


APPENDIX I 


SEA ICE FLOE SIZE DISTRIBUTION 
by 


D. A. Rothrock and A. S. Thorndike 


t 


Sea Ice Floe Size Distribution* 
by 

D. A. Rothrock and A. S. Thorndike 
Polar Science Center 
University of Washington 


Abstract 

Sea ice is broken into floes whose diameters range from meters to a 
hundred kilometers. This fragmentation aflects the resistance of the ice 
cover to deformation and the melting at floe sidewalls in summer. Floes are 
broken by waves and swell near the ice edge, and, throughout the pack, by 
isostatic imbalances, thermal cracking, winds and currents. In winter, they 
are welded together by freezing. 

Floe size can be measured by several properties p- for instance, area 
or mean caliper diameter. Two definitions of floe size distribution seem 
particularly useful: F(p), the fraction of area covered by floes no smaller 
than p\ and AT(p), the number of fioes per unit area no smaller thanp. A 
theorem from stereology states that F{p) can be measured by sets other 
than areas, such as the fraction of a line or of a point set covered by floes no 
smaller than p. If N behaves like p a for small p, where p is mean caliper 
diameter, a must be greater than -2 so that the small floes occupy finite 
area. If — 2<a<-l, the perimeter of small floes is infinite. 

Several summertime distributions have been measured. On a log-log 
graph, their slopes (local values of a) range from -1.7 to -2.5. One 
distribution follows a power law; the others have steeper slopes for larger 
floes, and more gradual slopes for smaller floes. Another sampling strategy 
is to measure the lengths of line segments on floes. The distribution of 
these chord lengths is equivalent to the distribution of floe diameters. 

The variance of an estimate of the fraction g of area covered by floes in 
any size range is p(l-p)A' -1 , where K is the equivalent number of 
independent samples. K can be found from the autocovariance of the 
indicator function for the chosen size range. For line sampling of a narrow 
range of floe diameters, K is the ratio of the sample length to the floe 
diameter. 


§1. Introduction 

The sea ice covering the Arctic Ocean is not a uniform continuous sheet like the ice 
that might cover a small lake. Instead it has irregular top and bottom surfaces and is 
broken into distinct pieces, called floes. In the summer these floes are-easily identified in 
remote images of the ice pack such as Figure 1; they Eire somewhat rounded in shape and 
are separated from each other by a lacey region of open water. They have diameters 
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ranging up to 100 kilometers. In the winter, floes still exist but they are less easily 
identified. Then the ice pack appears to be highly fractured, but the cracks, called leads, 
often cannot be resolved where the floes are actually in contact. To make matters worse, 
floes are constantly being "welded" together as new ice forms in the leads that separate 
them. 

The division of the ice pack into floes affects its large scale geophysical properties. It 
seems likely, for instance, that the way the ice pack deforms in response to forces applied 
by the winds and the currents and at the coastline is controlled more by the geometry of 
the pack than by the material properties of sea ice. In this paper we discuss some of the 
geometric properties of the ice pack including the sizes and shapes of the ice floes, the 
arrangement of floes in space, and the abundance of floes of various sizes. 

Physical processes that determine the floe geometry include: failure under 

horizontal forces applied by the winds and currents, failure during isostatic adjustment, 
thermal cracking, flexural failure in the presence of ocean swell, melting around the 
perimeter of floes, and abrasion with adjacent floes. These processes break large floes 
into little ones and cause floes to decrease in size. Buring the freezing season these 
processes are roughly balanced by the tendency of adjacent floes to be welded together. 
At present it is not known which processes are most important. Perhaps an 
understanding of the geometry of floes and how the geometry changes during the annual 
cycle will stimulate research on the governing physical processes. In any event, an 
understanding of the geometry of the ice pack is of interest in its own right for a number 
of practical applications associated with transportation in ice-covered seas and with the 
design of offshore structures intended to survive in the presence of ice. 

The first section of this paper defines the floe size distribution, clarifying the notions 
of "size" and "distribution". The second section describes techniques for measuring the 
floe size distribution. A sampling theory in Section 3 indicates how many measurements 
must be taken to resolve the floe size distribution to a specified accuracy. Some 
measurements that test the sampling theory and others that illustrate floe size 
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distributions of different ice covers are presented in Section 4. 

Measurements are presented primarily to illustrate points of technique or approach, 
and to indicate the range of distributions in different seasons. These have been made 
from three sources: 

1) an aerial photographic mosaic of the region of the Arctic Ice Dynamics Joint 
Experiment (AIDJEX), 18 August 1975 (Hall, 1979), reproduced here as Figure 1, 

2) part of a LANDSAT summer image, 185 km x 155 km, 18 August 1973 (No. E-1391- 
22283), 

3) four U-2 aerial photographs of areas about 30 km on a side, taken June 21, 1974 
(flight No. 74-101; frames 9635, 9637, 9639, and 9641). 

The only other published observations we are av are of are in Weeks et al. (1980). Those 
measurements were made from airborne X-band side-looking radar data. We also mention 
a project currently underway at the Scott Polar Research Institute to measure floe size 
distribution near the ice edge from satellite imagery (A. Cowan, personal 
communication). Other work on this topic is reported by Losev (1972). 

§2. Floe Size 

Let P{x,y,t) represent some measurable scalar property of the ice pack at the 
location (x,y) and at time t . In other contexts one might take P to be the ice thickness, 
the surface temperature, or some other local variable. But for our purposes we will 
consider P to be related to the size of the floe covering the point ( x,y ) at time t . If no 
floe covers the point, we take P = 0. Several measures of size are of interest: 

1) area: The area of the x,y plane covered by the floe in question. 

2) diameter of the largest inscribed circle: This is not as useful a measure of size 

as the total area of the floe, because it is not particularly easy to measure. 

3) mean caliper diameter: Imagine calipers consisting of two parallel lines. A 

caliper diameter is the distance between the lines when each touches one side of 
the floe without penetrating the interior. The average of these readings as the 
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ealipers are rotated through all orientations is the mean caliper diameter. 

4) perimeter: The length of the curve defining the floe boundary, assuming the 
curve is smooth and its length finite. This quantity may be related to melting on 
the sidewalls of floes. 

We now show that these properties are highly correlated so that a measurement of 
any one property gives an approximate value of any other. To study the relations between 
properties we selected the AIDJEX summer mosaic and digitized the perimeters of all floes 
with diameters over about one kilometer (Figure 2). All of the properties defined above 
were calculated for these 782 floes. The distribution of area/(mean caliper diameter) 2 , 
for example, is shown in Figure 3. Subsamples of 45 large floes and 379 small floes have 
substantially the same distribution. 

Table 1 summarizes these relationships; for example, a floe with a mean caliper 
diameter p has an area of 0.66 p 2 ± 0.05 p 2 . Note how closely the relationship 


Table 1. Statistics on certain ratios of floe properties 



Disc 


Floes (Summer Mosaic) 




Mean 

Standard cr/ Mean 

Deviation ftr) 

Area 

(Mean Caliper Diameter) 2 

•7- = 0.785 
4 

0.66 

0.C5 

0.08 

Inscribed Circle Diameter 
Mean Caliper Diameter 

1 

0.77 

0.09 

0.12 

Perimeter 

Mean Caliper Diameter 

7T = 3.14 

3.17 

0.04, 

0.01 

Area 

(Inscribed Circle Diameter) 2 

— = 0.785 
4 

1.14 

0.21 

0.18 

Perimeter 

(Inscribed Circle Diameter, 2 

7T = 3.14 

4.19 

0.54 

0.13 

Area 

(Perimeter) 2 

(47T)- 1 = 0.080 

0.065 

0.005 

0.08 
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perimeter = tt • mean caliper diameter 

is satisfied.' This expression is exact for any convex shape (Miles, 1978); a small concavity 
adds perimeter without altering the mean caliper diameter. The floes in this sample were 
in an intermediate stage of break-up. Floe shapes at the end of the melt season and in 
winter may be slightly different. 

§3. floe Size Distribution 

Within some geographic region S (e.g., the Beaufort Sea) at some instant t (e.g„ the 
first day of August), we consider the scalar property P to be a random function of space 

P(z,y;t) (x,y)eS ,t€T 

The observed function can be thought of as having been drawn at random from an 
ensemble consisting of all possible functions. 

Any particular realization of P contains a full description of the F-geometry. In 
general, this geometry is very complicated, and much of it is irrelevant since it would not 
be repeated in different realizations. Our aim is to define statistical properties of the F- 
geometry that are common among different realizations and therefore describe the 
ensemble from which the realizations are drawn. Many more realizations of ice floe 
geometry will have to be studied before this goal can be reached. The present emphasis 
is to develop the general statistical concepts, not to document the full range of possible 
geometries. 

We assume that the process is ergodic in whatever statistic is of interest to us. This 
means that one realization over a sufficiently large region S can teach us all we could 
learn by looking at numerous realizations. 

Fractional Area 

Suppose P is a property of an ice floe, so that P is constant over any particular floe. 
Then one can define the fractional area F in S covered by floes for which the property P 
is no less than p. Formally, use the Heavyside function H(q) = 1 if g § 0, and 0 if q < 0, 
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to write 


F(P,S) ~ ^JjHlPip,y)-p ]dxdy . (1) 

iS denotes both a region of the plane and its area without confusion. This definition avoids 
problems with "edge effects". If a floe with P = Po lies partially in only that part of its 
area within 5 1 is counted, but the value of p D pertains to the whole floe. F(p,S ) is a 
decreasing function of p. If p is a property that cannot be negative, then 


limJf(p,5) =Jf’(0 + ,iS) 
r-*o 

is the concentration of ice. F(0,S) equals unity. 


Thus far, 5 has been taken to be an area, but one might choose to sample floe sizes 
by some other set of lower dimension: either a line drawn through the image, or a set of 
regularly or randomly placed points. Sampling in such sub-spaces is called stereology 
and has countless applications. Early papers, for instance, concern estimating corpuscles 
in tissue and particulates in minerals (Wicksell, 1925; Nicholson, 1970). The subject 
contains a theorem that shows how measurements from different sample spaces are 
equivalent (Miles, 1978). The theorem states that if some subset U of S is selected, and 
U'cU is the subset of U for which P exceeds p , then the ratio of the measure of U’ to the 
measure of U is, on average, equal to F{p,S) 


where E denotes the expected value or ensemble average operator, and fi is the 
appropriate measure (the number of points, the length of line or the area). 


Number density 

A second way to describe the distribution of floe size is to find the number of floes 
per unit area of a region 5 for which P is no less than p 

N(p.S) . 

Again, edge effects are best handled by counting the fraction of a floe's area contained in 
S, while retaining the value P of the whole floe. N(p,S ) is a decreasing function of p. 


The total number of floes in 5 is 5 ‘N (0,5), which may be finite or infinite. 

The number of floes larger than p and the fractional area covered by floes larger 
than p are related by 

F(p,S) = ~fa(p)dN(p) 

p 

where a(p) is the area of a floe with property^? (taken, say, from Table 1), and -dN(p) is 
the number of floes with p<P<p+dp. Use of one rather than the other is a matter of 
choice. The advantages of F are its bounded behavior for small p and its appearance in 
the fundamental theorem of stereology. 

It might seem useful to consider the total number of floes larger than some cut-off 
po, N(po,S ) and to define the fraction of these floes larger than p 

s{p ' p ° )= fw • pipu ■ 

Q is a cumulative probability function which increases to unity as p decreases to p 0 . Its 
usefulness is limited, because the information about hew densely the floes are packed in 
space has been lost, and because Q depends on po as strongly as on p. If we envision 
picking floes at random, this probability function assigns each floe equal probability, 
whereas the quantity F, taken as a cumulative probability, weights each floe (and the 
open water) by its area. 

Power law distributions. 

Consider the distribution of mean caliper diameter p. For small p, N can behave like 
p a , ot^O with ot=0 corresponding to the case of a finite number of floes in 5 . For a<0, the 
number of floes smaller than any finite p D is infinite. Of course, the total area of such 
floes must remain finite 

Pa Pa 

- / p z dN = —a f p a+x dp < oo 
o o 

which is satisfied only for a>—2. The total perimeter of the floes is proportional to the 
integral of pdN, which is finite only when a>— 1. The three possible cases are: 


a = U, finite number of floes, 

— l<a<0, infinite number of floes, finite perimeter, 

— 2<a<— 1, infinite number of floes, infinite perimeter 

Over the range for which the power law holds, the floe geometry has a certain self 
similarity 

= g(Jc) . 

That is, for small floes, floes of a fixed size ratio occur in numbers of a fixed ratio. In this 
case there is no natural length scale to the geometry; provided that we choose a sub- 
scene which has small scale structure and does not fall entirely on a single floe, the small 
scale structure looks the same under arbitrary magnification. It is a common experience 
to confuse sea ice images with quite different scales, as in Figure 4, and that’s why. Self 
similarity is illustrated by the Apollonian gasket in Figure 5, which is characterized by a 
power law distribution N ~p -1,307 . 

Exponential distributions. 

A construction with quite different properties is shown in Figure 6. It is called a 
Poisson field, because the perpendicular distances from the (randomly oriented) lines to 
an arbitrary origin are distributed as a Poisson process (Solomon, 1978). Unlike the 
Apollonian gasket, it has a finite number of pieces in a finite region S. Circles inscribed in 
the pieces are distributed exponentially 

N(p) = E{Q)e~ p/x 

where p is circle diameter, and X is the length scale of the pieces. The Poisson field bears 
only a slight resemblance to winter ice, and none to summer ice; the infinite length of the 
construction lines precludes there being many small floes in the space surrounded by a 
few large floes. We do not expect to observe exponential distributions in nature. 
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Techniqu.es and observations. 

One technique we tried involved the use of a set of discs, graduated in size from 3 
mm up to 2.5 cm, corresponding to floes with diameters between about 1 and 10 km. 
Beginning with the largest disc, one can identify all the floes into which the disc fits. 
These floes have an inscribed circle diameter 5 no smaller than the disc diameter. 
Working with k discs, one can establish k points (N,5) of the number distribution of 
inscribed circle diameter. We used five discs. Floes more than half in the test area were 
counted. The method is appealing in its simplicity and is sufficiently accurate to 
distinguish the range Df distributions occurring in nature. 

This method was applied to the LANDSAT image and four U-2 photographs, the 
distributions for two of which are shown in Figure 7. Inscribed circle diameter was 
converted to mean caliper diameter p, using <5 - (0.77 )p from Table 1. 

Another technique is to digitize each floe boundary on an x,y digitizing tablet by 
tracing a cursor around the boundary of each floe. Some precaution must be taken to 
digitize all floes in the chosen size range and lying even partially within the test area. 
Each floe is represented by the coordinates of twenty to fifty points around its boundary, 
from which one can easily calculate all floe properties. 

This method allows one to treat edge effects without error, counting the (areal) 
fraction of a floe within the test area while measuring the floe property of the whole floe. 
The test area for the AJDJEX summer mosaic is the box shown in Figures 1 and 2. 
Digitizing requires equipment and some set-up time - for instance, writing programs to 
compute floe properties — but in an operational mode, this method is probably as quick as 
using disc cutouts and provides accurate measurements of all floe properties. 

The AIDJEX summer mosaic was measured hi this manner. Its number distribution is 
shown in Figure 7, and its fractional area distribution in Figure B. 

The mosaic and the LANDSAT image are data from the same geographic area in mid- 
August, two years apart. The difference between the two distributions (in Figure 7) is 
evidence of considerable interannual variability. 
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Striking spatial variability is also present. The distributions denoted by U2-35 and 
U2-39 in Figure ? were observed at the same time from regions only 50 km apart. 
Further, one can interpret the ice in U2-39 as being a broken version of the ice in 112-35. 
To the extent that this interpretation is valid, we see that only the large floes broke; U2-39 
has fewer 10 km floes but more 1 km floes. 


§4- Chord length distribution 

In winter, when it is difficult to identify floes and assign values of P, piece sizes can 

CL 

be characterized by the distribution of chord lengths. At each point A on a random line, R, 
assign a chord length C(x) of the longest segment containing x and lying on a single piece 
of ice. 


Either end of the chord touches a crack or lead. From R select a random test 
segment S of length A, count the number of chords within S for which C is no less than c , 
and divide by L to define the chord distribution M(c ). For a chord only partly in S, count 
the fraction of its length in S', but assign a value C equal to its total length. 


The chord and floe size distributions are closely related. In this instance, take the 
floe property to be the mean caliper diameter p. If a rectangular area with sides a and 6 
is sampled by a test line S parallel to side b , we have that 

, number of chords l r . probability that chord is between c and c+dc, l 
between c and c +dc ^ i ^ DB lers given that 5 hits a floe of diameter p I * 

( 2 ) 


, probability that S 
hits a floe with diameter p 

Using the definition 


number of floes with 
diameter between p and p+dp 


G(.c,p) - probability { chord > c IS hits a floe with diameter p $ 


m(c) = and n(p) = 

dc dp 


V 


and the densities 
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(2) becomes 


6 m{c)dc - f 


-#dc} 

£ 

dc j 

a 


(a 6 7 i{p)dp) 


or 


m(c)~ -f --^~ p -p n(p)dp . 


For circular floes with 


( 3 ) 


G = [l~4l 1/2 
P 

equation (3) is an Abel integral equation which can be inverted to give the diameter 
distribution in terms of the chord distribution 


M(p) = ~f (c 2 ~p 2 )~ l/ *dM(c) . (4) 

71 c=p 

Parts of this development follow Ripley (19B1, §9.4). 

Figure 9 shows measurements of both N(p) and M(c) from the AIDJEX summer 
mosaic, and a calculation of N{p) from M{c) observations using (4). Measurements of 
G(c,p) for real summer floes are shown in Figure 10. 


§5. Sampling the floe size distribution 

It may be inefficient or even impossible to determine the fractional area distribution 
F by measuring the area of each floe in a region R . In this section we discuss the errors 
that arise when F is estimated using a subset U off?. A general theory is presented and 
then applied to several specific sampling strategies where U is taken to be a box, a 
random straight line segment, or a regular lattice of points. 

Let P(x,y) be a measurable property at the point (z,y) and let F denote any 
condition on P. Define the indicator function II(r,7/) to be 1 if P satisfies T and zero 
otherwise. For example, if P represents mean caliper diameter, and T the condition 1 km 
< P < 2 lan, the indicator function will take the value II = 1 over all floes with diameters 
between one and two kilometers, and will be zero elsewhere. Use the symbol g to denote 


the fractional area of R for which P satisfies T. With our earlier point of view that P is a 
random function, we have 

g = f n(a:,3/) dR ~ Prob (P satisfies T) 

R 

Here the measure of R, fiR, is included to normalize g properly. The probability of g is 
assumed to be independent of x,y, any point in R has equal a priori probability of 
satisfying the condition F. 

Now let g be an estimate of g based on the random subset U . 

g = flu 1 £ n(x ,y) dU . (5) 

We now show that the mean of g is g and show how the variance of g depends on the size 
of the sample U. 

The mean and variance ofg. 

Since P is a random function, so is II. II has the trivial probability density; 
Prob(II=l) = g\ Prob(II=0) = 1 —g. Thus ETl equals g. Then using ( 5 ), Eg equals g . 

From ( 5 ), we have 

var g = E (g-g ) 2 = E^Ljj 1 f U{x ,y)dU -g 

= flu 2 Rn(x,y,x \y ') dUdU' , 

where /?n is the autocovariance function for the indicator function. Its properties are 
discussed in the next section. This equation describes how the variance of the estimator 
depends on the sample set U and on the way the floes satisfying the condition T are 
arranged in space, as characterized by Rjx. 

The autocovariance of the indicator function. 

The indicator function II depends on the property P being studied and on the 
condition F. In another context, P might be associated with ice thickness and F with any 
narrow range of thicknesses (Rothrock, in press). But here we associate P with mean 


fiv 1 £n(x\y')dU' - g 


( 6 ) 
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caliper diameter and T with any narrow range of diameters. In both cases, II has certain 
spatial structure that shows up in its auto covariance. 

Consider the one dimensional analogy in which the line is broken into intervals of 
random lengths, and the property P is interval length. Let T be the condition 
D < P(x) <D + e. Then II(x) is a random function of x taking the value zero for most x 
and the value 1 over intervals of length greater than D and less than D + e. The notable 
structure of II in this case is that all of its "up" segments have length D . For small e the 
up segments will be rare and the intervals between them will have random lengths. The 
situation is approximately that of a single "up" segment of length D placed at random on 
a line segment of length L»D . To determine the autocovariance for this random 
function, note that 


J?n(a) = £[IT(x+a) -n/L] • [ll(x) - D/L]=ETL{x+a)Ti(x)-D z /L z . 

The expression -ET^x+ajntx) represents the probability that x+a and x both lie on the 
random up segment that we denote by V. We have 


■fi’IIfc+ajnfx) = probability [(x+a) € V and x € F] 

= probability (x e V) - probability [(x+a) € V | x e V] 

This last probability is evaluated by considering all possible locations that allow V to cover 
x, and the fraction of those that also cover x+a. Then£ , n(x+a)II(x) is 


nrob ability (x € 7) ■ over * ap segments centered on x and onx+q 

p y ' * length of segment centered on x 

l-|a | /D , |a|<D , 

-9 ] 

0 , |a|>D . 


For large L , this is approximately i?n(°0, but we adjust the amplitude so that the variance 
/?n(0) is 9 )• keeping zero for large a. Then 


l-|aj /D , 


i? n (a)Kp(l-p) 


0 


\a\<D 

[a|>2J 
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In two dimensions, a similar argument can be used to find-AnC*.* 1 ) = i?n(s) where s is 
the distance between x and x\ Take the floes to be discs of diameter D , and again denote 
one rare disc by V. For V to cover x, its center must lie within a circle of diameter D , 
centered on x We have that Fn(x')n(x) is equal to the product of g and the fraction 

overlap of circles whose centers are a distance s apart 
area of circle centered on x 

which is 


2iT'h(§) 


where 


h(x) = 


ccs -1 ®— | x | (1— a; 2 ) 


0 


\x\<l , 

M>i • 


Again adjusting the variance, we have 


^n(s) = ff(l”ff)2Tr -1 /i(J-) . (7) 

figure (11) allows comparison of equation (7) with observed autocovariance functions for 
three different narrow ranges of floe size. 


Sampling by area 

Returning now to the question of the variance of g, consider the case where the 
sample set U is a box of side A. Then 

A A A A 

var g - A~*f f f f -R nte.l/.s'.y ') dxdydx'dy' 

0 0 0 0 

Make the coordinate change x',y' *+ s,6 with s = [(a;— a;') 2 + (j/-!/') 2 ] 1 ^ and 8 — tan -1 
{y'-y)/{x'-x). With s and 8 fixed, the region for integration for x and y is the overlap of 
two boxes of side A whose lower left-hand comers are offset a distance s in the 8 
direction: 

viT Bn -a | cosS | 4-s|sine| 

varg = A~* Iff J R?i(s)s ds dL8dx dy 
£=o <5=o *=o y=o 

'JZA 

— A~* J’ (27e 4 2 - BAs +2s 2 )sRn(s)cfe . (8) 




miTiT I.T 


To illustrate the use of this result, substitute the approximate expression (7) for 
into (8) to obtain 

'JEi 

varg = g (l-g)Z-n~ i A~ 4 f (27iA 2 -8As + 2s 8 )sA(^-)ds 

o " 

y/skD - 1 

~ g( t l~-g)27T~ 1 D z A' mZ f (Zn — BDA^ij + ZD z A~ s r] z )7]h(j])d7) 

o 

For A »D , this is roughly 

varg = g(l~g){~irD z )A~ 2 ' . (9) 

In many equations that follow, the variance of g is proportional to g (1—g): 

varg ~ g(l-g)K~ l , 

For M independent point samples (Bernoulli trials with probability of success equal to g), 
the constant of proportionality is simply K — M. For area sampling we have that 
K ^A z /{rD z /A) — (box areaj/fdise area) which we can think of as an equivalent number 
of independent samples. If, for example, we wish to estimate g for a class of floes that are 
thought to cover about 0.10 of the area, and we want the estimate (5) to be good to ±0.02, 
then varg = (0.02) a = (0.l)(Q.9)iT' 1 , which implies K = 225. This can be achieved by 
selecting 225 point measurements separated widely enough to be independent or by 
measuring densely over a box of side A = (n/A) l/z K l/z D , roughly if 1 ''' 2 or 15 times the 
diameter of the floes in question. 

Equation (8) web tested using the digitized AIDJEX summer mosaic. Using the 
measured autocovariance of the indicator function for a fairly narrow range of floe size, 
varg in (B) was evaluated for several values of A, as shown in the right-hand column of 
Table 2. An independent measure of the variance of g was obtained directly by choosing 
several different test areas of side A, measuring g in each one, and noting the mean and 
variance of the set of samples of g. These are shown in Table 2, and confirm the 
theoretical results. 
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Sampling along a line 

Let U be a random line segment of length L. At each point along the line, the floe 
diameter (not the chord length) is determined and g is evaluated using equation (5). 
Then the variance is 

L L L 

varg = L~ Z J CRji(\x—x'\)±c dx’ = 2L~ 2 f (L~s)Rji(s)ds 
oo o 

Using the form given above for i2n and taking L3>D , we get the approximation 

varg=g(l-g)4:n~ l f(l-x)h(xj-)dxBig(l~g)j- , L»D . (10) 

Sampling on a regular lattice 

When U is an MxM lattice of points with spacing A, the expression for the variance of 
g becomes 

var g = 

i=ys ii=u=i 

where S^i is the distance between the ij and kl lattice points. The covariances ) 

can be thought of an an M z xM s matrix of covariances between any pair of lattice points. 
Hie diagonal elements are simply g(l—g). If A is chosen large enough so that the points 
are uncorrelated, then the non-diagonal elements are zero and 

varg = g(l-g}M~ 2 . (ll) 

For a smaller lattice spacing some of the non-diagonal covariances will be positive and the 
variance of g will increase. 

To summarize, we estimate the fractional area covered by ice floes of a certain range 
of sizes using equation (5). The variance of this estimator can be found using the auto- 
covariance iJn of the indicator function in equation (6) for any sampling strategy. The 
variance is inversely proportional to K, the equivalent number of independent samples, 
which is given by the approximations 



i X 


i 1 


i 

i 


i 

i 

i 
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area sampling: 
line sampling: 
point sampling: 


K ra area of sample / area of floes being considered 
K « length of sample / diameter of floe 
K rs number of points separated by at least a floe diameter 


§6. Conclusions 

We have devoted this preliminary discussion of the floe size distribution of sea ice to 
questions of definition and of measurement, with emphasis on the errors which arise in 
estimating the distribution from a limited sample. Our preferred definition is the number 
distribution N(p) giving the number per unit area of floes whose mean caliper diameter 
exceeds p. N(p) increases as p decreases to zero. In some of the data sets we have 
examined, N behaves approximately like p° with -1.7 < a < -2.5, but we see no reason to 
expect a power few or any other simple analytical form to be valid for all p. We find 
changes in the distribution from year to year and from one region to another. 

The theory given for the sampling errors is summarized in equations (9), (10) and 
(11). The central idea is this: the error in an estimate of the abundance of floes of area a 
in a sample of area .4 depends on the number of independent samples, which is roughly 
A /a. 

The several floe properties related to floe size that we considered, mean caliper 
diameter, perimeter, area, chord length, diameter of inscribed circle, appear to be 
roughly equivalent. The choice of what property to measure may be made based on the 
kind of data and measuring systems available. For manual measurements, we 
recommend sampling the chord lengths along random lines. 

The important geophysical problems are to relate the geometric properties of the ice 
pack to its mechanical and thermal properties, and to relate changes in the geometry to 
the mechanical and thermal forcing. Floe size, ice thickness, and surface roughness are 
three geometric properties which have been studied. Others of potential interest are the 


ii 
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•widths, lengths and branching properties of leads. 
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Figure Captions 

Figure 1. Mosaic of aerial photographs of summer pack ice (Hall, 1979). The dark areas 
are open water between floes. The box measures 95 by 110 km. 

Figure 2. Digitized floe boundaries in the summer mosaic (Figure 1). 

Figure 3. The distribution of the ratios of area to squared mean caliper diameter. 

Figure 4. Views of sea ice on two different scales. The frame on the left is about 2B km 
wide, that on the right, 2.4 km. The texture on the left of the high level image 
is due to clouds. Melt ponds can be seen in the lower altitude photograph. 

Figure 5. An Apollonian gasket (from Mandelbrot, 1977). The disc sizes are distributed 
as a power law. 

Figure 6. A Poisson field, illustrating a construction with an exponential distribution of 
piece sizes. 

Figure 7. The cumulative number distribution of mean caliper diameter N(p). For the 
solid lines, mean caliper diameter was directly measured; for the dashed 
curves, it was calculated from measured inscribed circle diameters. 

Figure 8. The cumulative area distribution of mean caliper diameter F(p) for the 
summer mosaic. 

Figure 9. (above) The cumulative number distribution of mean caliper diameter N(p) 
for the summer mosaic: (•) measured directly, and (+) calculated from 
measured chord distribution (solid line below). The solid chord distribution 
was smoothed from the actual observations (*). 

Figure 10. The density of the ratio of chord to mean caliper diameter for the floes in the 
summer mosaic and for discs {-dD/dy, where y = c/p). 

Figure 11. The autocovariance of the indicator function. The solid line is the theoretical 
expression (7) for rare discs. The symbols show observations from the 
digitized summer mosaic for several narrow ranges of diameters. 
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° f aeria * photographs of summer pack 
The dark areas are open water between floes 
95 by 110 km. 
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Views of sea ice on two different scales. Frame 4a is 
about 28 km wide.. The texture on the left is due to 
clouds . 
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Figure 8, The cumulative area distribution of mean caliper diameter 
F(p) for the summer mosaic. 
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Figure 9. (above) The cumulative number distribution of 

mean caliper diameter N( p) for the summer mosaic: 
(•) measured directly, and (+) calculated from 
measured chord distribution (solid line below) . 

The solid chord distribution was smoothed from the 
actual observations (•). 
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Figure 10. The density of the ratio of chord to mean caliper 

diameter for the floes in the summer mosaic and for 
discs (- dG/Dy , where y “ c/p). 
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Figure 11. The autocovariance of the indicator function. The 

solid line is the theoretical expression (7) for rare 
discs. The symbols show observations from the digitized 
summer mosaic for several narrow ranges of diameters. 
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ICE THICKNESS DISTRIBUTION— MEASUREMENT AND THEORY 


by D. A. Rothrock 


INTRODUCTION 

In considering the historical literature relevant to the subject of ice 
thickness distribution, I was struck by the paper by Wittmann and Schule 
(1966) at the Lake Arrowhead symposium on the arctic heat budget. They 
presented a synthesis of visual observations from aircraft showing ice 
concentration, lead spacings, percent coverage of several ice types, and 
ridging indices in the subregions of the Arctic Ocean. Their format 
emphasizes the basic elements of present thinking — first, that we want 
properties averaged over space scales which include many ice features, and 
second, that these average properties vary on yet longer scales over the basin 
and from one season to the next, in ways of interest to the climatologist. 

Wittmann and Schule showed what data were available from operational 
observations just as Soviet scientists were attempting to forecast ice con- 
centration (Nikiforov et al., 1967; Doronin, 1970). In the intervening 
fifteen years, effort has been directed to defining what data we want for 
climatological and forecast modeling of sea ice, and why we want them. One 
data set ve want is ice thickness distribution. A major reason is that it 
tells us the mass of ice in the ice cover. Other reasons are the subject of 
other chapters. Suffice it to say that ice thickness and its distribution 
influence components of the surface heat balance, the salt balance of the 
upper ocean, light penetration and biological productivity, and traff inability 
on top of, through, and beneath the ice cover. It is unfortunate that 
thickness alone does not determine these phenomena; they are affected by snow 
cover, temperature, ice formation history and surface relief • 


Presented at the NATO Advanced Study Institute on Air-Sea-Ice Interaction, 
September 27-October 10, 1981, Maratea, Italy. 
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so that writing some flux as a function only of thickness sweeps a lot under 
the rug. 

THICKNESS AMD ITS DISTRIBUTION DEFINED 

^t-^w>3?sp ^ ir y p ficg 3r ■» £ gooph ygjj&sM&et ^1 though one can draw a slab of 
ice or a cross section of a ridge schematically with several straight lines, 
when one sets about measuring thickness, it is no longer such a clear-cut 
notion. There are cavities in ridged ice — some connected to the air or 
water. Ponds of fresh or brackish water covered with a skim of ice can lie 
in the troughs of both the upper and lower surfaces. How should a slush of 
ice crystals be treated? We suppose that one has satisfactory conventions 
for dealing with these phenomena and that there is a top surface of the ice 
(not the snow) at a height h^ measured up from the sea surface, and a bottom 
surface at draft h^ measured down from the sea surface. The thickness is 

h = i, t + V 

From thickness measurements which resolve ridges and leads, we want 
to know average properties over regions hundreds of kilometers in diameter. 

A 

For instance, to evaluate mean thickness H for a region R centered on (x,y), 
we write 


H(x,y) - x 




JJ 

eR 


h(x f ,y’)dx ? dy’ , (where A = dx'dy*) 


‘II 


as an integral, over thickness, of the area covered by each thickness 


II* 


da(h) _ 

A = 

R all 


h g(h) dh 




JS-- XXX-- —F 
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The ice thickness distribution g(h)dh is defined as the fraction of R with 
thickness between h and h + dh. If we want, say, an average heat flux () 
and Q depends on h, we can again integrate over all thicknesses 

Q = X jj QfM 2 -' >y t )]dx , dy f s | Q(h) g(h) dh 
R all h 


Since many fluxes at the air-sea-ice interface are strongly related to thick- 
ness, it appears that measurements of thickness distribution would have many 
applications. 

Two observed distributions of draft (which is nearly the same as thick- 
ness) are illustrated in Figure 1. The Beaufort Sea distribution shows 5% 
coverage of 0-to-l meter ice and 40% coverage of 2-to-3 meter ice. On the 
other side of the ocean, the Fram Strait ice is much heavier, with roughly 
five times as much thick ice as in the Beaufort Sea ice. The 2-to-3 meter 

category is comparatively empty with only 3.5%. Thu^.the Fram Strait ice 

vt’frU* flT it-J cffi-crruTa 

has a greater mass, is probab l y frictiona lly roughe^, and yet with its 10% 
of 0-to-l meter ice is producing ice faster than the ice in the Beaufort Sea 


sample. 

Although the distributions in Figure 1 look strikingly different, the 
question must be answered: are these real geophysical differences or could 

these two distributions be different samples from the ’’same" ice cover? An 
empirical answer is obtained by taking several samples from each area, and 
noting whether all samples from one region are similar to each other and 
different from the samples from another region. 
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An analytical technique for answering questions about sampling errors 
is rooted in probability theory. We adopt the notion that h£^, is a stochastic 
process; the thickness h at a fixed position ^gpis a random variable. In 
this framework, the thickness distribution is a probability function 


b 

p - | g(h)dh = probability that a<h<b Cl) 

a 

CO 

The expected value of thickness is the mean of the distribution H = / gdh. 

0 

By taking the process to be nearly stationary (varying only on scales of 
hundreds of kilometers) and ergodic, we can estimate statistics of the process, 
such as the mean, or distribution, or spectrum, from spatial averages of data. 
In addition, we can find confidence intervals for these estimates. 

MEASURING THICKNESS AND ESTIMATING ITS STATISTICS 

There are two aspects to measurements of thickness statistics: measuring 

thickness accurately at a point and sampling it at enough points to estimate 
a statistic satisfactorily. All methods are seriously limited either by 
measurement errors or sampling errors or both. 


Submarine Sonar 

The most satisfactory measurements of thickness distribution are obtained 
from upward-looking submarine sonar. Table 1 lists several cruises from 
which data have been published. The sonar does not directly measure thickness, 
but the distance from the transducer mounted on the submarine to the closest 
point of the bottom surface within the cone of the sonar beam. A pressure 
sensor measures the transducer depth. The difference between these distances 
is some approximation of the ice draft. 
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Table 1. Submarine sonar data. 

Beam width 

Submarine 

depth 



Submarine 

Period 

Location 

along track 

uncertainty 

Length 

Reference 

Nautilus, 
Skate 3 Sorgo, 
Sea.dra.gon 

1957-62 

whole 

Arctic 


2.0 m 

>10 4 km 

Lyon 

(1963,1966) 
LeShack (1980) 

Dreadnought 

March 71 

86°-90°N 

6°-7°E 

o 

o 

0.3 m 

1000km 

Williams et al. 
(1975) 

Sovereign 

Oct 76 

80°~90°N 

25 o E-70°W 

17° 

0.2 m 

4000km 

Wadhams (1981) 

Gurnard 

Apr 76 

70°-77°N 

138°-155°W 

3° 

0.1 m 

1400km 

Wadhams and 
Horne (1978) 


There are several sources of measurement error. First, the depth of the 
transducer is not precisely known; in the more modern data the accuracy is 
about ±0.2 m or 0.3 m. The error is not so critical for thick ice, but it 
might make an estimate of the coverage of thin ice useless. It would be useful 
to know more about the properties of this error. Wadhams (1981) takes it to 
be white noise with zero mean and variance (0.2 m) 2 . 

Another error is due to the non-linear smoothing of the surface by the 
sonar beam — always increasing the observed draft. The problem is akin to 
surveying the ice with a leveling rod 15 m in diameter. Because no information 
is sensed from some troughs, there is no way to reconstruct the actual draft, 
although an approximate reconstruction can be attempted (Williams et al. 9 
1975; Wadhrns , 1981). An alternative approach is to simulate numerically the 
smoothing of the wide beam (17°) sonar on a profile measured with a narrower 
beam. The wide beam causes the draft distribution to be shifted to the right. 
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underestimating thin ice and overestimating deep drafts and the mean draft 
by about 10% (Williams et al., 1975; Rothrock and Thorndike, 1980; Wadhams, 
1981). 

The desired quantity is thickness, not draft. So^ although isostatic 
balance is achieved only over tens of meters, it is assumed to hold at each 
point in order to convert h^fx) to h(x) by 

v*> 


The error thus introduced has zero mean and is probably not important. 

To estimate a statistic of thickness, we integrate over some length L 
of the profile. For instance 

L 

H = ^ | h(x)dx 
0 

is an estimate of the true mean thickness H. If L is chosen too small, the 

a 

estimate is poor. The estimator H is a random variable with a distribution. 
Its mean is H: it is an unbiased estimate. The important quantity Is the 

A A 

variance of H which tells how close H is likely to be to the true mean. The 
variance is 




L 


1 

ji 


A A 

var H = E(H - H) 2 = E - 


£ [ hCx^dxj - H 



► 

h(x 2 )dx 2 - B 

* 



0 


t 

) 

• 


= ~ 2 JJ E JhCxj) - hJ jh(x 2 ) - Hjdx 1 dx 2 

00 
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E is the expected value operator. The integrand is 


the autocovariance 


of h, C h ; assuming the thickness is homogeneous and isotropic, it is a func- 


tion of separation only^ s = Jxj - X2I 


The double integral reduces to 


var H = 


I MH 


(s)ds 


by transforming to new variables s and t = + x £ and integrating over t. 

For the Gwpnoord profile whose autocovariance is shown in figure 2, this sampling 
error is listed in Table 2. A typical mean thickness of 3-to~4 m can be 
sampled by a record ten- to- twenty kilometers long with a sampling error of 10%. 


Table 2. The sampling error of an estimate of mean thickness from a line 
record of length L. 


L, km 

(var H) 1/2 

10 

0.38 

50 

0.17 

100 

0.12 

500 

0.05 

1000 

'0.04 


Although no error model has been developed to treat these several measure- 
ment errors and the sampling error as a whole, it is likely that mean thickness 
estimates from 100 km sonar profiles have an uncertainty of 0.2 to 0.4 m. 

How large is the sampling error for distributions shown in Figure 2? In 
particular, how accurate is an estimate of p, for the bin a<h<b? To answer 
this question, we need to characterize the spatial arrangement of the regions 
contributing to this bin. These regions are identified by the indicator function 
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p 

1, a < h(x,t) < b 
0, else 
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The regions where it - 1 are the regions we must sample to estimate p. Given 
profile data (and so, ignoring the y.and t dependence), our estimate of p is 
* the fraction of the profile intersecting these regions 

L 

p “ j; J ir(x)dx 
0 


The estimate is unbiased, since 



Ett dx 


= probability 


L 

1 f 

[a<h(x)<b] • — 

0 


dx 


P 


assuming the process is spatially stationary. The remaining issue is: how 

large is the variance of p? We have 


var p = E(p - p) 2 

f L 




1 


= 17 


LL 


00 


C lr (xi,X2)dx 1 dx 2 
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where 


C Tr (xi,x 2 ) 



Is the autocovariance of the indicator function. Stationarity provides that 
depends only on the distance between points s = jxi - X 2 J, and the double 
integral reduces, as before, to 

L 

varp=| J(l -f]V s)ds (2 > 

0 

'^he autocovariance has magnitude p(l-p) at s = O' [since Eif(x)i:(x) = 

Ex(x) - p], and is zero at s = ® [since Eii(x)ir(x + s) + p 2 as s + »]. If 
we define the correlation coefficient 


r(s) 


C (s) 
pCI-p) 


(2) becomes 


var p 



( 3 ) 


where 


£ = 




r(s)ds 


In practice we are interested in the case where r(s) drops off in a distance 
much smaller than L; then 


£ “ 2 J r(s)ds 

0 
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< *“ m, > 5, we find » within M* of i„.) ft. guantity ^ is the ^ 

correlation length scale; it is also a measure of the characteristic length 
Of individual profile segments with v = 1. 

The first portion of Table 3 gives values of 1, eowputed from the Beaufort 
Sea profile. Values as small as 10 m have considerable uncertainty because 
the horizontal resolution of the data is several meters. The resulting 

sampling error for p is also shown there and in Figure 3. The length scale 
Tor thin ice is much greater than that for thick ice. because thin ice is spread 
out m leads and Pdynyas, whereas thick ice is clumped in ridges. In the 
table, we see that the thinnest ice and the 5-to-6 meter ice are egually plenti- 
ful, but that the error cr. is an order of magnitude smaller for the thick er 
-e bin. To get an egually good estimate of the thin ice would reguire a 
record 5600 km long (although not necessarily in a straight line). 



Bin fa) 


Beaufort Sea T. „ An ^ 

A 

P 


0 to 

1 to 

2 to 

3 to 

4 to 

5 to 

6 to 

7 to 

8 to 


1 

2 

3 

4 

5 

6 

7 

8 
9 


l 

CO 

(m) 

1300 

230 

170 

34 

20 

14 

11 

10 

12 


9 to 10 n 


.053 

.089 

.42 

.20 

.087 

.051 

.031 

.021 

.014 

.011 




.033 

.018 

.026 

.0095 

.0051 

.0034 

.0023 

.0019 

.0017 

.0014 


yp 

(%) 

62 

20 


Fram Stra it L = .100 

A 

P 


8 

9 

12 

13 


.105 

.026 

.036 

.137 

.110 

.084 

.073 

.073 

.073 

.057 


0 

P 


.035 

.0076 

.0077 

.0063 

.0044 

.0033 

.0027 

.0026 

.0028 

.0024 


cr-j/p 

P 


33 

29 

21 

5 

4 

4 

4 

4 

4 

4 
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Assuming £ is not very different for another mature ice cover, the 

CD 

values of SL a from the Beaufort Sea profile have been used in (3) to evaluate 
the sampling error for the Fram Strait distribution (in Table 3 and in Figure 
3). Within the sampling error, the percentage of thin ice (0-to-l meter) in 
he two distributions is marginally distinguishable, but in all other bins, 
the distributions are measureably different. 

Heat balance estimates are sensitive to the thinnest 5 or 10 centi- 
meters of ice«--£*ayjra& t t hin vol - u asej . Given the measurement errors and 
sampling difficulty, present submarine sonar records cannot provide good 
estimates of such narrow categories of thin ice. 

Boring Holes 

The most direct way to measure thickness is to bore a hole or cut a 
core. Thickness can be measured to whatever accuracy is appropriate to the 
condition of the surfaces. So there is negligible measurement error. 

It can take a few minutes’ work to bore a hole. A reasonable strategy 
for estimating thickness statistics might be to determine the thickness h.. 
by boring n holes. To be independent the samples must be roughly 1 km apart; 
this is the spacing at which the autocorrelation function of h^fx) has its 
first zero (Figure 2). How close is the estimate of mean thickness 


i=l 


to the true mean H? The variance of H is 


var H *= o 2 /n 


n 


(var 


23 ' I 

62 I 

560 I 

2240 I 
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separation should be several kilometers (several times if we are sampling 
the 0-to-l meter bin. The situation reduces then to n independent Bernoulli 
trials with chance of success p; the variance is simply p(l-p)/n. Table 5 
shows that to resolve the distribution into ten bins each with ten percent 
of the area to an accuracy of 30% (that is, p * 0.10 ± 0.03) would require 
one hundred holes. T his. is . getting impra c tica l-. Boring might be practical 

i 

in applications which require a fairly coarse resolution of the distribution 
with generous tolerances. 


Table 5. Number of independent samples n inquired to estimate fraction of 
area p with standard deviation o eoual to some fraction of p. 

p - 


p 

n for 

a = 0.3 p 

__P 

n for 

a = 0*. 
P 

0.3 

26 

233 

0.2 

44 

400 

0.1 

100 

900 

0.05 

211 

1900 

0.01 

1100 

9900 


Other Methods 

There is a class of observing techniques that may eventually provide 
measurements of ice thickness from aircraft or satellite. Some are discussed 
in another chapter, but not necessarily with regard to sensing ice thickness. 
All are in an early stage of development as thickness sensing instruments. 
Included are radio echo sounding, visual and infrared photography, various 
microwave and radar sensors, and laser altimeters. Only radio echo sounding 
senses both the upper and lower surface. Altimeters sense the height of the 
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snow surface; to estimate the ice surface h t and then multiply by 

ice density ] “ 1 O 0 
water density J 

to obtain thickness h introduces unsatisfactory errors. The microwave radio- 
meters sense radiative properties at or near the top surface — properties only 
weakly tied to ice thickness for thick ice. 

Some measurements of thin ice fraction have been obtained from LANDSAT 
visual images with accompanying ground truth. The area of gray ice is measured 
from an image; the thickness of the gray ice is determined by coring or drilling. 
This method has provided data for tests of thickness distribution theory 
(Thorndike, 1980) . 

THICKNESS DISTRIBUTION THEORY 

Having observed the thickness distribution and having at least a rough 
idea of the thermal and ridging processes which determine it, the ice physi- 
cist has naturally been tempted to formulate a mathematical theory to describe 
the distribution and its evolution. Such a theory was introduced as part of 
AIDJEX modeling work (Thorndike et al., 1975), and is gradually being intro- 
duced into forecasting and climatological models (e.g., Hibler, 1980). Its 
use does involve several difficulties. First, since the distribution is 
difficult to observe, often no initial conditions are available, and model 
results cannot be verified (or refuted) . Second, largely for the same 
reason, model parameters describing the ridging process are not well known, 
so that the model results can only be considered rough approximations. And 
third, generally one is modeling the ice cover as a function of two (horizontal) 
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space coordinates and of time, and the addition of another independent 
variable, ice thickness, considerably enlarges the computational burden. 
Being thus forewarned that the theory is no panacea, we will nonetheless 
review it and see what strengths we find in it. 

The equation for the time rate of change of thickness distribution 
g(h,x,t) is 


3t 


" div ^ ” 3h (fs) 


({) 


(4a) 


The four terms on the right hand side, discussed one at a time below, repre- 
sent opening and ridging, flux divergence, ablation and accretion at the upper 
or lower surface, and lateral melting. The Lagrangian form is 


- ifi - g div u - -^-(fg) + <f> 


(4b) 


These equations are primarily statements relating geometry and kinematics. 

The real physics of the ice cover enters in the structure of ifi, f P and <J>. All 
three of these functions depend on the thickness distribution itself, making 
(4) a non-linear functional equation. Early work with this equation treated 
only the feedback of the thickness distribution g through ridging and opening 
More recent work (Maykut, in preparation; Hibler, 1980) has included the 
feedback by which summer melting (4> and f) is determined by the amount of radia- 
tion absorbed and hence the open water. 

To consider flux divergence, we focus on an arbitrary region R and the 
exchange of h-ice (ice of thickness between h and h+dh) across its boundary 
cur-ve C. The flux of h-ice across a line with normal n is n . ug dh. The 


total flux into or out of R is written first as the integral of the flux 
around C 

r 

(i n * ug dh dl 
C 

and then, by Gauss’ theorem, as 

| div(ug dh)da 
R 

This term accounts for the net loss of h-ice from R 

ff (gdh)da 
dC J 

R 

Since R is arbitrary, and fixed in space, the equation takes the differential form 

l 

•|| = -div (ug) 

This flux divergence term has two parts: -g div u, denoting the influx of 

h-ice by convergence, and -u * grad g denoting advection. 

Thermal growth f = dh/dt carries ice to a new thickness, in just the 
same way as velocity u = dx/dt carries ice to a new position, giving a term 
completely analogous to the flux divergence term 



This is a mathematical statement. For a more physical view of thermal growth, 
it is best to turn to the cumulative distribution G(h, t) . The governing 
principal is that ice thinner than h at time t is thinner than h + fdt at 
time t + dt giving 

G(h,t) = G(h + fdt, t + dt) 

Expanding the right hand side about (h,t) gives 

G(h,t) = G(h,t) + • fdt + |“ dt + 0(dt 2 ) 

or 

j!G _ _3G 

3t 3h 

which is the integral form of (5) . 

Lateral melting if eats away at the edges of floes to reduce their area. 
The process conserves area: an amount of open water is formed equal to the 

total ice-covered area melted so 

+ = sJfiCh) - A(h)j 

where s depends on the available heat (but not on h) , 6 is the Dirac 

00 

delta function and &(h) is positive with J £(h)dh = 1. The volume loss 
of h-ice shJl(h) is proportional to the area of floe edges which we can 
take to be proportional to hg(h); hence, £(hY'°^k(h), and the integral 
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constraint determines that £(h) ~ g(h). The unknown s is determined by 
stating that 

CO 

Q = pL | h<j>dh 
0 

where Q is the rate at which heat becomes available for lateral melting 
(per unit area of R), P is the density of ice, and L is the latent heat of 
fusion. Solving for s and substituting gives 


- 5(h) - gCh) 


where H = / hgdh. 

0 

Ridging and opening processes are embodied in the redistribution func- 

f ' ~ '*v / fJ-J*- 

tion on which there are two constraints. The ita ec har aen ljprocesses can 
only rearrange existing ice, ridging thin ice to produce thick, but cannot 
alter the mean ice thickness^ so by assumption 


CO 

I 


h$dh - 0 


Furthermore, the formation of new area of open water less the loss of area 
by ridging must exactly accommodate the area of ice imported by convergence. 
Hence, ignoring the thermodynamics (setting f *= $ = 0) and integrating (4) 
over all h gives 


CQ 

| ifidh = div u 
0 


t IX '."I ** — 


V, 
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The definition of i|) is hypothetical, relying on assumptions and 
reasoning about idealized situations (Thorndike et al., 1975; Hibler, 1980); 
there are few hard facts about redistributions 

f^Open water formation is represented by a delta function at h = 0, 
and ridging by a function ]i(h) which is a sink of thin ice and a source of 
thick ice. The degree to which each of these processes occurs in redis- 
tribution depends on strain rate invariants. We denote the sum of the 

• m 

principal strain rates by e_ (= div u) and the difference by e TT . These 
can be expressed in polar form as a module's |e[ = (e^. 2 + e^ 2 )^ , which 

m * 

tells the quantity of the deformation, and an angle 0 - tan 
which tells the quality of the deformation. 

The combination of ridging and opening into ip is assumed to be propor- 
tional to strain rate in the form 

^ = |e| { a o ( 0 ) <500 + a r (6) u(h)J 

The coefficients a and a need to be observed in nature — the idea is that 
o r 

divergence (small 0) will favor opening (larger smaller o^) and conversely 

that convergence (0 getting up toward tt) will favor ridging. But even in 

pure shear (0 - %) the separation of floes in some locations and coining 
- *2 

together of floes elsewhere will insure that both opening and ridging occur 
s imui taneous ly . 

Just how the ridging function p varies with h depends on the present 
state of the thickness distribution (since we don*t want to ridge ice that 
isn*t there), and on two parameters G* and k. G* is the cumulative fraction 
of the thickest ice assumed to ridge. It is taken to lie in the range 0.05 
to 0.15, allowing ice up to — 




about one meter to ridge. The second parameter k is the ratio of the thick- 
ness of ridged ice to its original thickness. Of course, even the notion 
that this ratio is a constant is a strong contradiction of reality. The value 
originally assumed (k = 5) seems to have been too small, producing no ice 
thicker than 6 m and allowing the pack to deform too easily. (Ice strength 
is thought to be tied to ridging and to increase with the size of ridges being 
formed.) A value of k - 15 is better in both respects. 

An alternative formulation of redistribution has been advanced in which 
ridged ice is quantified not by its thickness distribution but in terms of 
mean ridge height and ridge intensity (Bugden, 1979). The motivation, of 
course, is that ridge statistics are readily available from aircraft visual 
observations and laser profiles, so a model predicting them would be testable. 
This approach has its drawbacks, though, and has not been hotly pursued. 

Ridged ice is not the sheet threaded with long triangular rubble piles seen 
in schematic diagrams. In reality old and new ridges overlap, making the 
identification of a pressure ridge a matter of each observer’s particular 
definition. So ridge statistics are not uniquely related to areas or volumes 
of ridged ice. Hence an equation involving areas of undeformed ice and 
ridged statistics for deformed ice is less clean than the formulation involving 
only ice areas. The more desirable path is to state the theory wholly in 
terms of thickness distribution and improve techniques for extracting it from 
data. 

The first theoretical solutions for thickness distribution were calcu- 
lated for a single Lagrangian point from (4b); that is, space was not a 
variable. We will examine solutions of that problem here. Deformation rate 


is a function of time, known from sets of buoy tracks, often lacking frequen- 
cies greater than 1 cycle per day. Growth rate is a function of thickness 
and time, known from climatology, and lateral melting is neglected. Starting 
with an observed distribution as an initial condition, solutions change in 
time but hover near a realistic distribution. Thus the theory is compatible 
with our understanding of ice dynamics. 

The theory, of course, can do more than predict the distribution. It 
shows us which terms are important to the balance in various ranges of thick- 
ness and just how large these terms are. Table 6 shows the annual average 
of net thermal production (growth into the range less growth out) , ridging 
and open water gain, and import by convergence. The thinnest category is 
maintained by newly formed open water growing thermally into the next cate- 
gory; some ice is ridged out of this category. For the rest of the thin 
ice (0.1 to 1.6 tn) , the net thermal production is nearly balanced by ridging 
loss. For equilibrium ice (1.6 to 6.4 m) both thermodynamics and ridging 
provide sources. A balance is only made possible by the divergence. Thick 
ice (>6.4 m) is not in balance. It is gaining 0.7% per year because the 
ridging gain cannot be compensated for by this divergence and the net melt. 

Whereas ice thinner than the mean adjusts fairly rapidly to a given 
deformation and growth regime, it is probably futile to think about an equili- 
brium distribution for thick ice. For one thing, the determining variables 
are poorly known: the long term net divergence, the infrequent large con- 

vergences capable of making very thick ice, and the melt rate and rate of 
erosion for thick ice. For another, the adjustment time for thick ice is 
slow: it takes thirty years for twenty-five meter ice to approach equilibrium 

thickness. Only the Beaufort Sea may contain ice that old. 




Table 6. An average over one year of the terms in the thickness distribution 
equation (4b) . Each term is also integrated over a thickness 
range. Units of terms are percent per day, of thickness distri- 
bution, percent. 


Terms in (4b) 


0-0.1 m 

0.1-1. 6 m 

X. 6-6*4 m 

6.4-25.6 m 

Thermal flux 
gradient 

ah (fs) 

-0.59 

0.43 

0.013 

-0.003 

Open water 
production 

A 5 

0.73 




Ridging loss 

- C *L 

-0.14 

-0.44 

0 

0 

Ridging gain 

C 

0.00 

0.01 

0.015 

0.014 

Import 

- g div u 

0.00 

0.00 

-0.028 

-0.009 

Rate of change 
of g 


0.00 

0.00 

0.000 

0.002 

Thickness 

distribution 

b 

f 

gdh 

j 

2.9 

20.0 

58.9 

18.1 


a 


The balance in Table 6 is no surprise but simply shows that the model 
behaves much as the ice it was designed to mimic: open water grows to thin 

Cf c-1. 

ice and is ridged into thick ice which continually ablates^ The theory eluci- 
dates some unforeseen details, and most importantly tells the rates of these 
various transfers and shows how nearly balanced they are. 


i 
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A WORD ON NUMERICAL TECHNIQUES 

One numerical procedure for solving this problem involves integrating the 
functional differential equation 

A # A 

H - j *dh - G div u, where « £ + f 
0 Dt 

for each time step along characteristics satisfying 
f (h,t) 

and then to interpolate the new G(h + fdt, t + dt) to some fixed h-grid 
G(h, t + dt). A large area of open water freezing at the same time (late 
summer appears as a step in G(h) or a spike in g(h). This step cannot be 
resolved by any fixed grid. But to make matters worse, repeated interpola- 
tion at each time step continues to smooth out the step. The numerical solu- 
tion then is only a shadow of the correct solution. An alternative procedure 
(Colony, personal communication) is to allow the h-grid to float, moving with 
the characteristics. A step can be followed by two nearly identical charac- 
teristics on which the values of G are quite distinct. Figure 4 shows solu- 

*s 

tions using Colony s characteristic grid, a coarse fixed grid (a la Thorndike 
et al. , 1975) , and a much denser fixed grid (153 h-values at spacings varying 
from 2 to 50 cm). Table 7 shows further comparisons. The integration starts 
with all thin ice, proceeding in time steps of one day for a year; the defor- 
mation rates are those computed from the AIDJEX manned camps (Colony, 1978). 
The values k = 15 and G* - 0.1 were used, and A and C were those identified 
with the 30° teardrop in Rothrock and Hall (1975). 
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Table 7. Variables computed by different numerical procedures. 


Numerical 

procedure 

mean 

thickness 

thermal 

production 

coverage 

by 

ice < 1 m 

coverage by 
ice between 
1 and 2 m 

coverage 

by 

ice > 6a 


h, m 

fdt, m 

G(l), % 

G(2)-G(l) ,% 

1-G(6) , % 

Characteristic grid 

144 

163 

23 

69 

95 

Dense fixed grid 

144 

171 

34 

58 

95 

Coarse fixed grid 

176 

246 

45 

23 

96 


The large step at 1.2 a is the remains of the initial thin ice. The 
smaller step at 18 cm is from the freeze up several days before the integration 
stops. These are clearly not going to be captured by a coarse fixed grid. It 
is more surprising that the finer fixed grid smears the large step out over 
nearly one meter. Another curiosity is that even the coarse fixed grid repre- 
sents the thick ice better than the thin ice from which it is produced. I 
see two lessons here. First, special numerical care is warranted for these 
discontinuous functions. Second, some variables are robust against numerical 
errors, and some are quite sensitive; it would seem prudent to know into which 
class one’s favorite variable falls. 

G can also have steps in space. If there is ice of thickness H to the 
left of an ice edge at x = X, G has a jump along the curve shown in Figure 5. 
Representing this function G(h,x) accurately will require more than single 
values on a fixed (h,x)-grid. And as H and X change with time, the jump in 
G will be lost if special care is not taken. 


TESTING THE THEORY 


Several tests of the theory (Nye, 1975; Thorndike et al., 1975; Rothrock 
and Hall, 1975; Coon et al., 1977; Thorndike, 1980) have focused primarily on 
the dynamics of thin ice and open water and have compared these quantities to 
satellite observations. They show that: 

i) the theory adequately simulates thin ice with the chosen forms 

of A, C, and p , and the value G* 0.05 to 0.10, 

Li 

ii) the considerable uncertainty in deformations estimated from the 
motion of only several paints can cause considerable discrepancy 
between computations and observations, in particular of summer 
open water, 

iii) the value k = 5 produces ridged ice that is insufficiently thick. 

As another test of the theory, focussing particularly on the thick ice, 
the following situation has been simulated. Ice forms in the Siberian marginal 
ice zone, and drifts across the ocean. It passes the pole in roughly three-to- 
five years, judging from the drift tracks of North Pole 6 and 10 and the Sedov, 
Starting on September first with an initial condition of 100% open water, the 
thickness distribution model just discussed (k = 15, G* = 0.1, 30° teardrop, 
AIDJEX manned camp deformations, climatological growth rates, no lateral 
melting) was integrated for several years. 

Some results are compared in Table 8 with submarine data from the vicinity 
of the pole. The theory has not produced enough thick ice. That was a short- 
coming in the calculation by Thorndike et al. (1975), but there k was 5, and 
no ice over 6 m was produced. Here, k is 15, and ice is produced in the whole 
thick end of the distribution, up to 23 m. Furthermore, the shape of the thick 
end of the distribution is similar to observations. Ice of the right thickness 
is being made but in too small an amount. The more likely cause for the dis- 
crepancy is in the deformations driving the ridging. They were observed in a 





ature ice cover in the Beaufort Sea, and may be quite inappropriate for young 
ice in the marginal ice zone. 


Table 8. Comparison of theory to data in a transpolar drift. 



Mean 

thickness 

Coverage of 
ice thinner 
than 2.2 m 

Coverage of 
ice between 
2.2 and 5.6 m 

Coverage of 
ice thicker 
than 5.6 m 


h, m 

% 

% 

% 

Theory, 
after 3 years, 
October 

3.0 

74 

12 

14 

Theory, 
after 5 years, 
October 

3.5 

64 

16 

20 

Sovereign, 
sec. 17-21 
October 

4.4 

13 

59 

33 

Dreadnought , 
March 

4.6 

12 

55 

33 


Sovereign and Dreadnought data are from Wadhams (1981) . Drafts are 
multiplied by 1.11 to obtain thickness. 


The open water at the start of the simulation has, after five years, 
grown to 1.9 m and covers 30% of the area. The Dreadnought thickness distri- 
bution shows a strong peak at 2.9 meters' thickness (Williams et al., 1975, 
Figure 10). The Sovereign peaks occur at slightly greater thicknesses. This 
deficiency of the model could be due to poor growth rates or to an underesti- 
mation of rafting and ridging of the young Siberian ice. 

The procedure of combining inputs and thickness distribution observations 
from different times and places is not going to provide any further improvement 


in the theory. As the ditty has it, "the'v gone about as fer as the* kin 
go." A more controlled and complete experiment is necessary, in which a 
Lagrangian region is monitored repeatedly by submarine and satellite, and 
the deformation history and thermal environment is documented continuously. 

WHITHER THICKNESS DISTRIBUTION? 

Several avenues are available for advancing our knowledge of ice thick- 
ness distribution. 

1) The parameters in the model are pinned down as well as available data 
allow. A controlled mass balance experiment would be valuable in which 
the model inputs and outputs were simultaneously monitored. 

2) Sampling theory can determine how large a sample is required to resolve 
the differences we think we observe between two distributions. This should 
be a prerequisite for a mass balance experiment. 

3) Remote sensing is unlikely to resolve the thickness of thick ice, but a 
program combining satellite estimates of thin ice and submarine sampling 
of thick ice would be valuable. 

4) Substantial concentration of a single thickness creates numerical problems 
that need more attention. 

5) The ice cover moves in pieces and not as a continuum as presently assumed 
in theory and in the analysis of deformation. A more direct approach to 
redistribution could be developed in terms of the movement of pieces. 
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FIGURE CAPTIONS 


Figure 1. Distributions of ice draft for two regions. The bin size is 

one meter. The Beaufort Sea data are from Gurnard profiles (the 60 km 
offshore profile of Rothrock and Thorndike, 1980), and the Fram Strait 
data are from Sovereign profiles (a 97 km profile, section no. 2 in 
Wadhams , 1980). 

Figure 2. Autocovariance of the Beaufort offshore profile (from Rothrock 
and Thorndike, 1980). 

* 

Figure 3. The range of probable values for the distributions in Figure 2. 

In each bin, the range p - Op to p + Op has been filled in — stippled for 
the Beaufort Sea distribution and crosshatched for Fram Strait. 

Figure 4. Cumulative thickness distribution computed from theory using three 
numerical schemes: characteristic grid, coarse fixed grid, and fine 

fixed grid. 

Figure 5. A line of discontinuity in G on the (x,h) plane. Open water lies 
to the right of x * X. 




one meter. The Beaufort Sea data are from Gurnard profiles 
(the 60 km offshore profile of Rothrock and Thorndike, 1980), 
and the Fram Strait data are from Sovereign profiles (a 97 km 
profile, section no. 2 in Wadhams, 1980). 
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Figure 2. Autocovariance of the Beaufort offshore profile (from Rothrock 
and Thorndike, 1980). 



VRfitT ( m) 

Figure 3. The range of probable values for the distributions in Figure 2. 

In each bin, the range p - (T p to p + r p has been filled in — stippled 
for the Beaufort Sea distribution and cross hatched for Fratn Strait. 
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Figure 5. A line of discontinuity in 
G on the (x f h) plane. Open water lies 
to She right of sj * X. 
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Kinematics is that branch of physical science dealing with motion 
itself and not with the causes of the motion. Attention is focused, in 
this chapter, on how sea ice moves. The deeper question of why sea ice 
moves as it does is addressed in the chapter on sea ice dynamics. 

For some scientific and practical applications, knowledge of the ice kine- 
matics (as opposed to the dynamics) is sufficient. As examples, consider 
gouging of the sea floor in shallow water, the transport of natural debris 
or pollutants, the impact loads on structures, navigability in icy regions, 
and the movement of the ice edge. Aside from these applications, the study 
of kinematics is a basic step toward understanding the dynamics. 

Figure 1 illustrates the idea that the motion of sea ice is the response 
of the ice pack to external forces. To understand the response of the ice 
one needs to monitor both the driving forces and the motion. Since we will 
discuss only the motion here, there will be the underlying ambiguity of 
whether the results inform about properties of the ice itself or about the 
driving forces. We will see for example that the ice motion is nearly non- 
divergent. But we should not conclude on this basis alone that the ice is 
an incompressible material. The correct explanation could be that the net 
driving force is itself nearly non-divergent. Because of this ambiguity 
in the interpretation of kinematic data, it is useful to think of the study 
of the kinematics as part and parcel of the broader study of sea ice dynamics. 

The fundamental kinematic notion is that an pieces of ice have an 
identity which is preserved in time. If we identify a piece of ice by 
making a small mark on the ice surface at position X at time o , the assumption 
is that at some later time *fe , the piece of ice as identified by that mark will 
be at a position X . This defines the position function x(i,X) with x(o,X) “X 
We do not mean to imply by this that ice floes retain their identity indefi- 
nitely, only that it is possible to track individual points. Two points 
originating on the same ice floe may very well wander apart in time, but we 
assume we could keep track of them. 

There are difficulties with the notion of a position function. It is 
natural to consider the space domain to be the two dimensional surface of 
the Arctic Ocean, say .21 • At *tr=o there will generally be some points in X 
which are not covered by ice, so is not defined for some X^X. 
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Further, ice is always being created and destroyed thermodynamically, so 
the position function x(^,X) is only defined for an interval of time. 

Finally, the idea that the small mark retains its identity is question- 
able. This is a common problem in the definition of the motion of a 
turbulent fluid. After a time the marked fluid has mixed with the unmarked 
and it is no longer possible to assign a position to the mark. Something 
of the same sort happens with sea ice. Despite these shortcomings, the 
notion of a position function is the best we have. Furthermore this notion 
corresponds exactly with most of the observations we have. We will use 
the term trajectory for the function where X is held constant. The 

initial position X serves to label ice particles and we refer to the 
trajectory of the particle X. 

5.1 Observations 

The most common observations of ice motion are of trajectories of ice 
particles consisting of a sequence of measurements (& L , L i) , <■ ' ‘ 1 ' W where 

5- l = + Ci 


The measurement error is represented by Ei . 

In passing we mention that other kinds of measurements of ice motion 
have been made. Hunkins (1967) and McPhee (1978) for instance inferred the 
ice velocity using current meters suspended from the ice into what they 
assumed to be essentially a static ocean. These data were used to study 
motion on the time scales of hours. An attempt to study motions on much 
shorter time scales using accelerometers was made by Craig (1972). 

Many techniques have been used to measure ice trajectories. These 
are summarized in Table I. Each technique has its good and bad features. 
Over the past decade most data have been obtained by satellite positioning 
which works somewhat as follows. Suppose a stable frequency •/ is trans- 
mitted by a device on the ice. The signal received at the rapidly moving 
satellite will have a frequency •/ which has been shifted by the Doppler 
effect: ■£■ = ^ + A-f . The Doppler shift // is related to the rate of 

change of distance between the device on the ice and the satellite 1 

j~H — #icet • The received frequency -f is measured at several 

times. If the satellite coordinates are known at these times, then each 
measurement produces an equation with unknowns and . Generally 

several measurements are made during the 10 to 20 minutes it takes for the 
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TABLE I 


Techniques for measuring sea ice motion 


technique 

basic measurement 

sampling 

rate 

accuracy 

comments 

Hunkins etal, 1971 
celestial navigation 

azimuth and elevation 
of heavenly todies 

1 per day 

10 3 m 

weather dependent 

Martin efc al,. 1978 
satellite navigation 

Doppler shift of stable 
transmitted signal 

30 per day 

30-500 m 


Thorndike, 1973 
acoustic tracking 

travel time from ice to 
fixed reference on the 
ocean floor 

up to 1 per 
minute 

5 m 

limited range, elaborate 
instrumentation, rela- 
tive position only 

Tucker et al., 198 O 
Tahata et al., 1980 
radar 

travel time from source on 
shore to target on ice 

no limit 

< 1 m 

limited to near shore 

Hall,. 1980 
Landsat imagery 

location of identifiable 
ice features 

erratic 

80 m 

weather dependent, good 
space resolution, 
poor time resolution 

Hall & Rothrock 
synthetic aperture 
radar imagery 
(SEASAT) 

location of identifiable 
ice features 

1 per 3 days 

80 m 

no system available at 
present, properties of 
future systems unclear 


satellite to pass by. The several equations are solved simultaneously for 
yt and %-y. . Precise positioning requires a stable transmit frequency, 
precise measurement of the received frequency and precise knowledge of the 
satellite coordinates. With care, errors can be controlled to the order of 
tens of meters, as in the best uses of the Wavy Navigation Satellite System. 

Satellite systems which serve primarily to relay data from automatic 
data platforms (or buoys) to data processing centers also determine the 
location of the transmitting data platform, using the same Doppler positioning 
principle. The ARGOS system currently on the NOAA-B satellite can relay 
data from up to 30 sensors and determine platform locations to an accuracy 
of a few hundred meters about ten times per day. Fully automatic platforms 
cost in the neighborhood of $6,000 with additional costs depending on the 
desired sensors. 

Satellite imaging systems can also be used to measure ice motions, pro- 
vided that some features on the ice can be identified in a sequence of images. 
Because of their all season and all weather capability and good resolution, 
imaging radar systems will probably be best. The basic angular resolution 
of these systems is approximately the ratio of the wavelength of the radar 
signal to the diameter of the antenna. To achieve an angular resolution of 
10~ 5 (10 meters at a range of 1,000 kilometers) with a wavelength of 25 cm 
(l. 2 GHz) requires an antenna 25 km in diameter. Although such large antennas 
cannot he constructed in space, it is possible to synthesize large antennas 
by using data from several points along the satellite's orbit. The determi- 
nations of the geographical position of an ice feature with the SEftSAT data 
contained errors of up to 3 kilometers (see Hall and Rothrock, 1981). Unless 
these errors can be reduced the data are not particularly valuable for 
measuring the displacements over intervals of a few days or less. However, 
the errors are highly correlated in space and are essentially eliminated in 
estimates of the spatial variability of the ice motion. 

The attractive feature of imaging radar systems is their potential to 
sample densely in space. Hall and Rothrock’ s work suggests that it will be 
possible to track roughly one feature per square kilometer, which will 
resolve most of the spatial structure of the field of motion. Techniques 
for extracting data from the images or from the raw data are still rather 
primitive. No doubt satisfactory automated techniques for identifying and 
tracking features will he developed when the need arises. At present there 
is no imaging radar system in space. Planning is underway for a system to 
he in operation perhaps by 1985. 
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5.1.1. Sources of data 

The earliest ice motion data are the trajectories of ships beset in 
the ice. These are followed by the trajectories of numerous Soviet and U.S. 
drifting research stations, and more recently by the trajectories of auto- 
matic data buoys. The following list of data sources is by no means complete. 


Experiment 

Reference 

Description 

T-3 

Hunkins et al., 1971 

Monthly positions May 1962- 
December 1970 

Beset ships and 
drifting stations 
before 1970 

Hastings, 1971 

A chart showing most prior 
trajectories, one point per 
month; includes : Fram, Jeanette, 
Maud, Sadko, Sedov, Tegetthoff, 
British Trans-arctic Exp., 
Alpha, Arlis I, Arlis II, 
Charlie, North Pole 1-20, T-3. 

AIDJEX 1972 

AIDJEX Staff, 1972 

100 km triangle, Beaufort Sea, 
March and April 1972. 

AIDJEX 

1975-76 

Thorndike and Cheung, 1977 

Manned camp positions and 
velocities tabulated at 6 hr 
intervals, daily buoy positions. 

1976-77 buoys 

Thorndike and Cheung, 1977 

Tabulated daily positions of 
Beaufort Sea buoys. 

1979 buoys 

Thorndike and Colony, 1980 

Tabulated daily positions, 

25 buoys, analyses of surface 
pressure plotted daily. 

1980 buoys 

Thorndike and Colony, 1981 

As above. 

Fram I, 1979 

Hunkins et al., 1979 

March-May 1979; 8L°N, 9°W 

LOEEX 79 

Popelar et al . , 1981 

Three stations, 100 km spacing 
April, May 1979; 88°-90°N. 


5.2 The general circulation 

The main features of the long term circulation. Figure 2, are the clock- 
wise circulation in the Beaufort Sea — and the motion of ice from the Siberian 
coasts across the North Pole and through the Greenland-Spitsbergen passage. 
Time honored nomenclature for these features are the Beaufort Gyre and the 
Transpolar Drift Stream. Some handy numbers for these long term features are: 

center of Beaufort Gyre: 80°N, 155°W, half way between Pt. Barrow, 

Alaska and the North Pole 


| 


V 

1 
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time to make 1 circuit: 5 years 

time to traverse Transpolar Drift Stream: 3 years 

area flux through Gr e enland-Spi tsber gen Passage: 

300 km 2 /day or about 20/5 of the area of the basin per year. 

This pattern of motion exists only as an average over several years. 

On shorter time scales there are departures from the long term pattern. Con- 
sider for example the trajectories plotted in figure 3 for the year 1979 and 
figure k for 1980. The trajectories are characteristically meandering and 
convoluted showing that on monthly time scales the ice motion differs markedly 
from the long term mean. On shorter time scales, not resolved in these figures, 
the motion is even more irregular. Notice the major anomaly which occurred in 
the summer of I960 when the motion of several buoys for several months was in 
the opposite sense from the long term mean. These departures of the actual 
motion from the long term pattern in most cases represent the response of the 
ice to the passage of atmospheric systems. 

If we ignore the forces and examine only the ice motion, the departures 
from the mean circulation appear aperiodic and chaotic. The departures can 
be thought of as random but they are not without structure. Our objective 
now is to clarify this structure. 

5.3 Ice velocity 

The ice velocity field can be defined by the relationship 

tr*0 ft 

Here the particle label ){ plays no role. 

This definition is meaningful only if the limit exists. From a practical 
point of view it is useful only if the limit exists and is approached when 
the interval ft decreases to the time interval T between observations. 

In figure 5 several sets of observations of ice motion are plotted, 
showing the variation of one coordinate of position versus time. Successive 
data sets divide the sampling interval t by 15 and improve space resolution 
by the same factor, giving a sequence of closer and closer perspectives on 
the motion. By constructing velocity estimates, U-jJ-i) * x(i,X)]/kz 

for k decreasing to 1, one can examine the limiting process. In the first 
two figures of the sequence, continues to change appreciably for small k . 

In the last two figures, h-fct-fc) at most times "fc becomes almost independent 

of k for small k , implying that on these time scales^ X ^ 15 minutes, 
the ice indeed possesses a velocity. 
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to a definite limit implies that 
the increment is proportional to t for small t . 

For random processes, a useful condition for convergence is mean squared 
convergence in which the variance of the increment becomes proportional 
to T 2 for small 1 . This can be readily tested for ice trajectories by 
plotting £(xiir%X) versus z on a log-log plot, as is done in 

figure i . The process is differentiable if the graph has slope of 2 for 
small t . 

From a practical point of view the process ceases to be differentiable 

at measurement intervals t for which the slope of the graph departs appre- 
► 

ciably from 2. The evidence in figures S and 6 Implies that the ice does 
have a velocity and that it can be resolved with a sampling interval of 
about three hours. 

Time averaged velocities can be defined without reference to the limit- 
ing process. Let **- cvw be the time averaged velocity at £ . 



This quantity u.t't.T'*/ and its properties depend on the duration T of the 
time averaging. For example, the variance of will in general 

be less then the variance of U.(*fc,Jv because the ~T average has suppressed 
contributions to the variance on shorter time scales. Also, from a dynamical 
point of view, the equations which satisfies should involve ~f as 

a parameter. Different physical processes may be responsible for determining 
u.(£»T *) for different values of 7~. 

Typical ice velocities range from 0 to 20 cm sec** 1 . An extreme velocity 
of 1^0 cm sec -1 has been observed. Two histograms of ice speed are shown in 
figure 7 , one corresponding to a full year of observations, the other 

restricted to summer observations. The winter data contained several periods 
of essentially zero motion. During the summer the ice was never observed to 
stop. 

In the following pages particular importance is attached to the time 
and space variability of the ice motion. 'We will often refer to the variance 
of velocity, £j r fu.-o) 1 ^(v-V^ ) = f\ This quantity has "been evaluated from many 
observations; it varies appreciably with season and with location. Table JT 
gives estimates of the velocity variance for each of the buoys shown in figure 3 . 


Convergence of jVW+T/O - 54 /x 
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TABLE II 


Variance of Velocity 
1979 data, 1 point per day 


Buoy ID 

velocity variance 

number of data points 

1901 

59 cm 2 sec** 2 

342 

1902 

48 

295 

1903 

50 

284 

1905 

210 

231 

1906 

81 

299 

1907 

37 

229 

1908 

174 

6l 

1909 

38 

267 

1911 

19 

62 

1913 

84 

31 4 

1914 

83 

282 

1915 

172 

31 

1916 

57 

172 

1918 

35 

283 

1917 

29 

219 

1920 

54 

293 

1923 

63 

302 

1924 

167 

72 

1925 

67 

280 

1926 

168 

59 

1927 

62 (median) 

311 • 
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5.3.1 Velocity time series POOR QUALITY 

The ice velocities deduced from trajectories should he interpreted 
from the Lagrangian point of view since the measurements are made following 
the material particle. If many trajectories are measured simultaneously 
it is possible to obtain by interpolation velocity time series at fixed 
Eulerian points. Although the Lagrangian description is more directly 
measurable, it has the drawback of sampling both the time and space varia- 
tions in velocity. The differences between the Lagrangian and Eulerian 
points of view is not of great practical importance in the motion of sea 
ice. This is because most of the variability in the ice motion is driven 
by the wind, and the space patterns in the wind field move across the basin 
so much faster than the ice moves, that all points, Lagrangian or Eulerian^ 
experience essentially identical stochastic forcing. Consequently, Lagrangian 
and Eulerian time series of sea ice velocities look just about the same. 

A typical velocity time series, taken from an ice station trajectory, 
say, will have a mean velocity usually less than 3 cm sec -1 , and it may show 
a trend. "When these effects are removed from the time series what remains 
are the fluctuations with time scales longer than the sampling interval 
and shorter than the length of record T . The AIDJEX position data for 
instance, serve to resolve fluctuations on time scales from a few hours 
up to a few months, ( t * 1 hour, T & 1 year). Over this range of time 
scales the ice velocity has a power spectral density as sketched in figure 8 . 

The power spectral density is plotted for positive and negative frequen- 
cies corresponding to counter-clockwise and clockwise rotations of the 
velocity vector. At the end of this section an algorithm is given for cal- 
culating the spectra of vector time series. We interpret the velocity 
vector time series as a complex time series with the Fourier 

decomposition 


1(^*0 s J 


tusk j 

dl(cc) a Jlo 


r* 

The power spectral density is the real function oCt^) - Q.Cui) a. (u>) , defined 

in the frequency range < to < -jfc. . 

The integral of the spectral density over this frequency range is the 
St 

variance a , and the integral over any frequency band is the part 
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of the total variance contributed by fluctuations with frequencies in that 
band. A useful way to summarize the information contained in the spectrum 
is to state the fraction of the total variance coming from frequencies 
greater than a certain value. 


Table III Fraction of total variance exceeding certain frequencies 


frequency: 

period: 

1 cycle per month 
1 month 

2 

2 weeks 

k 

1 week 

15 

2 days 

30 

1 day 

6o 

12 hours 


58% 

k5% 

3h% 

12# 

1% 

3% 



The equation of motion for sea ice balances the ice acceleration against 
air stress, water stress, pressure gradient forces due to the sloping sea 
surface, and internal ice stress gradients. While it is not the purpose 
here to examine the ice dynamics, it is useful to relate features of the 
ice velocity spectrum to these driving forces. In the central basin about 
15 % of the variance of the ice velocity can be explained by the local geo- 
strophic wind (Thorndike and Colony, 198£). In fact, the ice velocity fluctua- 
tions are roughly proportional to the local wind fluctuations. This implies 
that the spectra of the ice velocity and of the wind should have anproximately 
the same shapes (see figure ? ). 

The water stress dc-pends on the difference in velocity between the ice 
and the upper ocean. If the ocean is at rest the water stress is simply a 
drag opposing the ice velocity and %ts only effect is to reduce the ice' 
response near the inertial frequency. If the ocean is in motion, the ice 
will be carried along with it, and should acquire spectral traits similar 
to those of the ocean. Unfortunately the spectral signature of motion in 
the Arctic Ocean is not known. The long term circulation of the upper 
ocean appears to be similar to that of the ice, with a clockwise circula- 
tion in the Beaufort Sea and a transpolar current flowing from Siberia 
through the Greenland-Spitsbergen passage. It probably is not productive 
to ask whether the ice drives the long term ocean circulation or vice versa 
because, in the long term, the ice should be thought of as part of the upper 
ocean. In any case the long term ocean behavior does not affect the spectrum 
we are considering since we have subtracted out the long term mean velocity. 
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Fluctuations in ocean currents on shorter time scales have been observed 
but it is not yet possible to say how great their effect on the ice velocity 
may be. Some evidence for these variable currents is listed here. 

1. Monthly variations are apparent in the currents deduced from wiqd and 
ice motion observations by Thorndike and Colony, 1982-. 

2. Subsurface eddies. Manley* s (1981 ) study of subsurface eddies during 
AIDJEX described numerous features with velocity of about 50 cm sec- 1 , a 
length scale of 10 kilometers, and a depth of 50-200 meters. He found no 
expression of these structures in the ice motion. 

3. Inertial oscillations . The balance between the ice acceleration and 

the Coriolis force leads to oscillatory ice motion with a period which decreases 
from 12.77 hours at 70°N to 12 hours at 90°N. Inertial motions are always 
clockwise — hence the negative value for the frequency — in the northern Hemis- 
phere because the Coriolis force always accelerates the velocity to the right. 
Inertial motions in sea ice were first described by Hunkins (1967 ) and have 
received subsequent study by McPhee (1979) and Colony and Thorndike (1980). 

The amplitudes can reach 0.20 m sec 1 " 1 during summer when the ice pack is 
comparatively loose. Their effect on the summer ice velocity spectrum is 
striking (see figured). 

4. Tides . Tidal currents in the central basin are small because of the 
great depth. Theoretical estimates are in the range of 1-2 cm sec -1 . Over 
the shallow continental shelves the amplitudes are predicted to be at least 
an order of magnitude greater (Kowalik and Untersteiner, 197? ). These 
theoretical estimates are for the lunar semi-diurnal tide. Evidence from 
tide gauges around the basin summarized by Sverdrup (1926) implies that the 
lunar semi-diurnal (period 12.47 hr) and the solar semi-diurnal (12 hr) tides 
are the dominant tidal constituents in the Arctic Ocean. 

The prediction of a large amplitude of the tidal current over the shelf 
is confirmed by Sverdrup's observations from the Maud ; Hansen may 

also have observed tidal motion in the ice pack surrounding the Pram. Because 
the tidal and inertial periods are so close it may be difficult to diagnose 
observed motions correctly. Nevertheless there are several differences 
between the two kinds of motions which can sometimes he used to distinguish 
them. First the tidal vector traces out an ellipse during one period. An 
ellipse can be viewed as the sum of a clockwise circle and counter-clockwise 
circle. Thus, unless the tidal ellipse should happen to be exactly a 
clockwise circle, it should have some expression on the counter-clockwise 
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side of the spectrum. Inertial motion is strictly clockwise. Second, at 
a fixed point, tidal motions should have a fixed phase. Inertial motion 
on the other hand, acquires a new phase every time the ice receives a 
sudden impulse of momentum. Third, the phase of tidal motion should vary 
smoothly and slowly in space. Inertial motions at different points may 
have no fixed phase relationships. 

The motion described by Sverdrup can hardly be mistaken for inertial 
motion. The smooth variations in phase as the Maud moved from Wrangel Island to 
the New Siberian Islands, and the clearly elliptical cycle described by the 
measured velocity vector are not compatible with inertial motions. 

Most ice trajectories which have been analyzed in the western literature 
have been over the deep basin where the tidal motion is small. Careful obser- 
vations and analysis would be required to detect a tidal component of order 
1 cm sec" 1 since it would be mixed with an inertial component which is often 
much larger and with a rich spectrum of other types of motions. 

The Soviet literature contains many references to tidal motion in sea 
ice. Doronin and Kheisin (197?) and Zubov (19^3) each devote several pages 
to the subject. 

A consequence of the tidal motion is the associated cycle of con- 
vergence and divergence caused by the difference in phase of the tidal 
cycle at different points. Periodic opening and closing of the ice in the 
shallow seas has an effect on the heat exchange between the atmosphere and 
the ocean and on the rate of ice production. The theoretical calculations 
of Kowalik and Untersteiner indicate maximum divergence rates in the shallow 
seas exceeding 10“ 6 sec _1 which is enough to produce one percent opening 
during the tidal cycle. Their theoretical estimates of the divergence rate 
associated with tides over the deep ocean are three orders of magnitude smaller. 

Calculation of cross power spectral density for two two-dimensional vector 
time series. 


Given: two discrete complex time series U-' L , V;, , t =V with sampling 

interval V . 

Step 1. Select M , the number of lagsj a good rule is N/5 . 

Step 2. Remove' the mean and .trend from each time series. 


Step 3. Compute cross covariances 

M-Jf 

i I v* Y 


Z u-up Vj* 

p-JL J 


X * A, , 


V. 


"5 = 


j- 

N 
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Step 4. Define spectral window 




\ ( J + ** 


?rJl \ 

M+t ) 


J?- - Mj - M 


Step 5« Calculate spectrum S and frequency cj 

-'irJLk&M 

S, - * L T * e 

k *«-M 


- Mj ' ■ ') M . 
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5 * 3-2 Velocity correlations in tine . 

An alternative characterization of a random function is its autocorre- 
lation function. The discussion in the previous section of the power spectral 
density of the time function at a fixed or moving point could have been 

given in terms of the autocorrelation function 

E >ftn) *■"(*> . JL E uli ^)u.ni) . 

£ art) uf-(t) 

The two functions Sfa) and RC't) are Fourier transforms of each other and 
therefore contain equivalent information. Tfhich description is the.;jmore useful 
depends on the application. The spectrum is useful for distinguishing physical 
processes with distinct characteristic frequencies; for example separating 
the free inertial oscillations from the wind forced motion. On the other 
hand the autocorrelation function is more useful for questions related to 
prediction or experiment design. For example the question: ,T How well can 
tomorrow's ice velocity he predicted on the basis of today's velocity,” has 
an answer involving the autocorrelation function Rtir) evaluated at T = 1 day. 

The autocorrelation function is complex. Its real part contains informa- 
tion about the lagged correlations of the fC component of velocity with itself 
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and the £ component with itself. The imaginary part has information about 
the lagged correlations between the SC and velocity components. If the 
velocity of a piece of ice undergoing an inertial oscillation has a positive 
pL component at a certain time, it will have a negative component three 
hours later. The imaginary part of the correlation should have negative 
peaks at 3 hours, 15 hours, etc. 

The autocorrelation function for the velocity time series at the AIDJEX 
ice station Caribou is Bhown in figure // . The autocorrelation function 
estimated by the 1979 buoy data for the central basin is shown in figure . 
The real part of the correlation falls to about 0.7 after one day, 0.4 after 
two days, and decreases slowly at longer lags. The large correlation at lag 
one day indicates that persistence (the forecast strategy which predicts 
that the future will be the same as the present) will have some success for 
one day forecasts. 

As expected the autocorrelation from Caribou has negative peaks in the 
imaginary part at 3 hours, 15 hours, etc. Generally though the imaginary 
part is small. If the inertial motions are not of interest in a particular 
application, the imaginary part of the autocorrelation function can safely 
be ignored. This is equivalent to treating the two velocity components as 
independent time series. 

5. 3. 2.1 Application of the time autocorrelation function . 

It is desired to estimate the ice velocity u. at time given observa- 
tions 2j, at times = . Suppose the observations have zero 

mean random errors £(■= -uCij) which have covariance 

. 

The delta function expresses the independence of errors at different times. 
Finally suppose the errors are independent of the actual velocities. Choose 
an estimator which is a linear sura of the observations. 

a 

It is desired to find the complex constants which give the best estimate 

It in the sense that E is minimized. This is 

accomplished by differentiating f with respect to the real and imaginary 


iV. ^ A.- - 
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parts of each °<j. and equating to zero. The result is the system of 
complex equations for the unknowns «,,•••.> o'-n : 

x 
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A further reduction is achieved by noting that 


E ?fl s E«j M * + 

Use of the matrix notation i 

gives the compact expressions 

£ = * 
M*A =P 


H SjW- " E u^u.* 

{jfj] , A' [*j] , M*{Eitju£ + 

BxfcJ 



which involve the autocorrelation function R and the measurement error 
variance . The estimation error is 





With v-hese two expressions we can answer a number of practical problems. 

1. What sampling rate is required to insure good interpolation to the 
mid point between two observations? For given sampling interval X and measure- 
ment error variance , we write eq. 1 as 



The solutions for the , which in this case are complex conjugates for 
reasons of symmetry, and for F are sketched in figure J3 . The solution 
technique extends trivially to interpolation involving more points . 

The optimal choice of the weights <*, and is not «*, - - } A . This 

linear interpolation between the two data points is the most natural scheme, 
and it is quite good for small X , and small o* 1 /^ , but it is not the 
optimal interpolation. Especially at large X , a smaller error variance 
is achieved by giving less weight to the observations. For very large X 
the optimal estimate is simply u = o since s in this case the 

two observations are so removed in time as to have no correlation with 
the desired velocity. 






13 


ORIGINAL PAGE tS 
OF POOR QUALITY 

1 

The measurement error O' affects and F slightly. The dependence 
of F on the sampling interval T is quite strong. It appears that' if 
one needs to know each velocity component to ± 3 cm sec”’ , F - F £u * = 

18 cm sec” * , a sampling interval of about six hours is required even with 
perfect measurements. 

2. What errors are expected in forecasts "based only on the present velocity? 


Here we obtain 


c< = 


__&{* y 

/ + crV$* 




ry _ , _ RtoKW 


The forecast error is plotted for two values of initial error in figure l4. 

The prediction error grows with time, approaching the velocity variance 
a?- for large 7 . Prom the sketch we see that an optimal 3 hour forecast 
captures about 75$ of the velocity variance, an optimal one day forecast 
captures about 35$ of the variance. 

The velocity was assumed to have zero mean in this analysis and^in 
constructing the autocorrelation function ( figures the long term mean 
was first removed from the data. In an actual application to a region where 
the mean is known the best forecast would be 

a • ■ u + - (2 ■ u ) • 


Note that the best forecast is not simple persistence ( ®< = l). For long 

(c&o) 

forecasts, the best estimate is the mean^ } as equation 1 correctly indicates. 
5. 3.2.1 Acceleration 

The time autocorrelation function for velocity, figures 11 and 12, behaves 
like h%/c for small X • It follows from this that the variance of the incre- 
ment in velocity is proportional to (not X?" ). This implies 

that the velocity is not time differentiable in the mean square sense. Attempts 
to measure accelerations precisely by sampling over shorter and shorter time 
intervals may merely result in larger and larger estimates of acceleration. Of 
course, time averaged accelerations exist and have finite variance. 


0^ = V«ur 


- 2fk(i-*vn) 

T is 2f/cT with e«/d&c. 


lU. 


which for small 


In earlier sections the emphasis has "been on the structure of the complex 
function U. of the real argument t , at a fixed or moving point. In this 
section, the emphasis is on U. as a function of position for fixed time. 
Certain results are more easily expressed if U- and pi are regarded as two 
dimensional vectors rather than as complex numbers. The autocorrelation 
function 




Cu.(/,)vu; E 


is a 2 x 2 matrix involving possibly four different functions, whereas the 
time autocorrelation function involves ^ust two. 

The structure of the spatial autocorrelation matrix for sea ice veloci- 
ties reveals some of the properties of the velocity field. To establish 
the basic ideas we first determine the structure of the autocorrelation matrix 
for a homogeneous, isotropic, two dimensional velocity field, borrowing from 
Kolmogroff’s (19^-1) classic discussion of three dimensional turbulence. 

A random field is said to be homogeneous if its statistical properties 
are invariant under a translation of coordinates, and isotropic if they are 
invariant under rotations and reflections of coordinates. In particular the 
autocorrelation matrix for a homegeneons isotropic field must satisfy 

R.(*i , .) ~ R(tf, . 7-t +• G. ) for any d , and 

for any orthonormal matrix M . Multiplication by an orthonormal matrix 
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accomplishes an arbitrary rotation or reflection of coordinates* Note 
that the change of coordinates is applied'to both the position and the 
velocity vectors. These definitions make precise the notions that the 
statistics do not depend on position or direction. By choosing the trans- 
lation vector IX or the rotation and reflection matrix suitably, we can 
see how these definitions of homogeneity and isotropy constrain the form of 
the autocorrelation matrix. 

First, fl( -X. depends only on the vector difference . To 

see this, take Cl = -/£/ . Note that this does not imply that R. depends only 

A* 

on the distance j* 

Second, R. is symmetric. This requires using a. * , M * "X and 

cl - %■> * to obtain this sequence of equalities: 

~ rJ- 

Ett- (njvfrj - E U.lo)v('^-pC i ) = j:(- IA (-0 )J(- v(x ,-■?!*)) - f tiCKjvCfi,) . 

Third, E"u.(a t o} v(% } o) = a . Here use the reflection invariance in the 
definition of isotropy by taking 



which changes the sign of the ij. axis. Then 

Euto, o)v(fto) = . Eb-fOj-ojl-vf^-o)] - - £ U(t> j0 ) V(# j0 J 
which is only true when £ uco,e) vC#,*) -O. 

Consequently the autocorrelation function evaluated at ( x > °) 

has the form 



involving only two function of distance. Bjj (-r) is the correlation between 
the components of velocity parallel to the line joining two points separated 
by a distance r. $]_(->-,) is the correlation between the components of 
velocity perpendicular to that line. 

More generally if X and %z are arbitrary points with f 2 ' + ($r 

06 s (JWGVV and s«X.B , then 
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” U 


A -- 


c&b -it- wS 

st^fe 


Then ) s y* - ^- z E A U.U) ^ * /4£>(r)A T c 

^ m tt>?8 Gi,(r)4-^eGL(r) ct^S ^6 (^r) - ^ \ 

\ mB a-6 (Ell W - GiCr)) a^.*s G]j W + §i-ry 


<0 


which is the general form for the autocorrelation matrix for a homogeneous, 
isotropic two dimensional velocity field. It involves two functions, j3j) 
and 13]^ which we will estimate below from observations. Further it involves 
the separation -t and the orientation 9 of the two points A and at 
which the autocorrelation is being evaluated. It may seem strange that 8 
enters the statistics for an isotropic field. The condition of isotropy 
does not imply that Fu.lo,o)U.i o) equals fu.(b,o) u.C D t" i ') however. A 
rotation of coordinates affects both the velocity and the position vectors. 

A valid statement is £ U.(o,o)u.(r f a) = £v(o,o) v(o, r) See figure 15. 

There is another constraint on the spatial correlation function. Let 
U.- t ~ (ut, V L ) be the velocities at any set of points %- L i and let 

Q be arbitrary real numbers. Then the linear sum 

s- £ Wu + % n 

i-/ 

has a variance which can be expressed in terms of ^ . The added constraint 
is that Vo-r (s)^o for all choices of } ef^ , ^ . This requires that 

the matrix consisting of the correlations' of the velocity component at any- 
set of points [^.| must be positive definite (all positive eigenvalues), and 
is equivalent to the condition that the two dimensional spatial power spectrum 
of velocity be positive for all pairs of wave numbers . 

Contradictions can arise if correlation functions are used which do not 
satisfy this constraint. For example, suppose we wish to consider a velocity 
field for which J3 f j (e) - l, 0]jfL) - 0-1 tM-itl Gj|(2L) - °*5 . The 

attempt to evaluate the variance of 5= dfq,&) - /.J u.(L,o) + a.yu.(2L j0 J i n terms 
of these correlations produces a negative result. Since variance is intrinsi- 
cally positive, we must conclude that no velocity field could have the presumed 
correlations. 
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f&P ' 



I 


As an e .ample of an analytical form for ^ which satisfies the' positive 
definite constraint and hears some resemblance to reality, define 


fiUf) r e ^ e 


- r/A y 


J\ ^QO)aw\ 


This example was constructed by taking a stream function jE with the posi- 


_ T JAt 

tive definite correlation function f \(x) z e 
spectrum is the Hankel transform of e 
which is always positive. Now define velocity components u. = -"Si/aji 
V = 


The corresponding power 
This is irrf ^ + ) / 


and deduce R.(r) - - 

w/||l » r 


and G, (r) - - . 

-i "3r l 


5 r 3. 3.1 Estimates of the autocorrelation functions. 

In order to estimate the autocorrelation matrix, simultaneous velocity 
observations are required at pairs of points separated by distances ranging 
from roughly 1 to 10 kilometers. We present data from the 1972 buoy pro- 
gram and from AIDJEX which define the functions and I3j_ only for 
distances exceeding 100 kilometers. The 1979 data were first interpolated 
to give velocity estimates at a uniform grid in space and time. 

tlyl -- U. . &* *■*$ = **> *» * 6 = 4 *' = t 

The mean velocity over all points i,j,k was removed. Then for lags 
the lagged correlations were found using 


ft 


JV_*f 

JVLIvLMt ^ ~ i-+*j J&j k 

3 i = i r> k-i 




For these calculations = 7, = k 3 and A£ = 200. The results are sketched 

in figure JC for the lags Jl ~ 1, ..., 7, w - M =0, and J? = n = 0, 

#7 = 1,..., L. The results show that F u.(cs,o; o) and E Vfo^oJ VLo,r) 

are similar functions of t* , end that F vCo,o) d-CrjO) , Fu.(c,d) f 

Ir u(b,o)v(r 1 o) j a.*.J E\/(t> J b)u.(0j7-J) are all smell. '.These observations 

are roughly consistent with the assumption that the field is homogeneous 
and isotropic. We will proceed then on the assumption, supported by these 
data, that the ice has a homogeneous, isotropic velocity field with the 
functions By and empirically determined from figure /4 . 

Correlations for the AIDJEX data were constructed by choosing two of 
the ice camps, resolving their velocities into components parallel and perpen- 
dicular to the line joining the camps, and correlating these components. 
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These correlations are also plotted in figure l£. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


5. 3*3. 2 Applications of the space correlation. 

The spatial autocorrelation function is the basis for sound experiment 
design and data processing. For instance, suppose observations of velocity u.l 
are made at /V points ;£ -with independent zero mean errors having vari- 

ance cr* . Estimates of velocity are to be made at arbitrary points pt. using 
these observations. What is the maximum estimation error? We use the tools 
developed earlier. To estimate the u. component of velocity at , we use 
the observations u-,j"-,4 h and construct the matrices 




where Rw.-VJ = ^ ^ j) etc. as given by Eq. ^ . The esti- 

mation error variance is given by Eq. Z . Figure // shows how the maximum 
estimation error depends on the separation between measurement stations and 
on the measurement accuracy. The curve labelled a*/^ = o refers to the idealized 
condition in which the measurement errors are zero. Thus it represents the 
interpolation error which is due to the intrinsic spatial variability of the 
velocity field. 


The standard deviation of the ice velocity itself is about ^*“10 cm sec -/ . 
A reasonable observational goal is to keep the interpolation errors below 2 or 
3 cm see -/ . This can be achieved if the raw measurements of velocity are 
good to about 3 cm sec -/ and the grid spacing is about LOO km. At larger 
grid spaeings the interpolation accuracy deteriorates rapidly. 

5.h Deformation 

The differences in velocity from place to place are responsible for the 
characteristic morphology of the ice. When the velocity difference between 
neighboring pieces of ice is such that they tend to move apart, a lead forms 
and widens exposing the ocean surface to the atmosphere. During the winter, 
ice growth is rapid over open leads. If the motion changes—perhaps because 
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of a change in the winds — so that these pieces of ice move toward each other, 
the open lead will close, and any new ice which formed there will need to 
rearrange itself to occupy a smaller area. Typical mechanisms for this 
rearrangement are rafting where one part of the new ice sheet overrides the 
other, and ridging where the ice is crushed into pieces which pile into 
ridges rising a meter or two above and sometimes many meters below the sur- 
rounding ice. Ridging and rafting are not restricted to-thinj.ee. If the 
closing motion continues, the original pieces of ice come in contact and 
one may override the other or, by grinding together, pieces may break off 
and pile up and down to form a ridge. The essential ideas here are that 
the ice accommodates divergent motion by increasing the area of open water 
rather than by stretching and thinning. It accommodates convergent motion 
by reducing the area of open water and by ridging and rafting. These processes 
link the morphology of the ice, characterized by such features as floes, 
ridges, and leads, to the spatial differences in the ice velocity. Ify - 
studying the spatial variations in velocity we may be able to understand 
better why the ice pack has the form it does. 

There is a second reason for studying the deformation of sea ice. As 
the ice pack deforms, stresses develop within it which tend to oppose the 
deformation. These stresses figure into the local balance of forces and 
therefore affect the motion of the ice. To be more precise, the balance 
of forces equation contains terns of the form where 

is the stress tensor. The stresses may in turn be related to the ice motion 
by a constitutive law 


stress = p (deformation) . 

One of the motivations for AIDJEX was to investigate the function f from 
a theoretical point of view and by using observations of the deformation and 
indirect estimates of the stress (see Rothrock et al. , 1980 ). 

The concepts which have been applied to the study of the spatial varia- 
tions in ice velocity are those appropriate to the analysis of the deformation 
of a continuum. The fundamental concepts are the partial derivatives of 

and with respect to pi and ^ . The line of thought is 

that a description of the large scale deformation in terms of large scale 
average derivatives should give some idea of how much opening and closing 
is going on on smaller scales and some idea of the state of stress. For 
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example under a purely divergent motion only leads are formed. During a 
purely convergent motion no leads are formed. On a small scale, shearing 
along a crack is of no morphological consequence. But large scale shearing 
is often expressed locally by opening at some cracks and closing at others. 

Mathematical formulation of these ideas requires two steps. Required 
first is an appropriate quantitative description of the large scale defor- 
mation in terms of observable parameters. The second requirements is for 
a functional relationship giving the opening, closing, or state of stress 
in terms of the observable parameters. 

Current practice is to describe the motion in terms of the large scale 
strain rate invariants, which themselves are constructed from the large scale 
velocity derivatives. The first invariant 


(T = — + 
Ft -w 




^b 


quantifies the large scale divergence and convergence 3 and the second invariant 

= fe -s) * b* + S/ J 


(lfl ,e)i where 

r& , *E*J* 


ifi 


expresses the total rate of deformation and 


0 -- 




0 $ & < V 


indicates whether the motion is predominantly divergence &~o , shear 
or convergence & '-‘W . 

In the theory presented by Thorndike et al. , 1975, opening and closing 
are assumed to be known functions of 9 and proportional to }£} . 


A = total opening = 

C = total closing = J Fj 


The functions ©<„ and e< r which give the total opening and closing in 
terms of the strain rate invariants may have a form somewhat as sketched in 
figure 18. Some constraints on the shapes of these functions have been 
deduced by Nye, 1975 and by Rothrock and Hall, 1975- From pp-.r? of Landsat 
images of the ice these investigators were able to measure the strain invariants 
and the total opening and closing. 
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The^f un cti onal relationships for A and C in terms of |F| and & 
can only he regarded as approximate since they do not take into account the 
actual geometry of the floes in the region being considered. Different 
regions with different geometries will respond differently to similar large 
scale strain rates. 

In practical applications of these ideas the procedure has been to 
take measured velocities at several points in a region of interest — perhaps 
the three corners of a 100 km triangle. Velocity derivatives are estimated 
from the measured velocities. The estimated velocity derivatives are com- 
bined to form the strain rate invariants JF) and 0 , and these are used 
to find the rates of opening and closing A and C , and the state of stress <T 

In the following sections we reexamine some of these ide*vs, particularly 
the notion of velocity derivatives, the influence of the ice pack geometry 
on opening and closing, and the difficulties in basing inferences about 
these quantities on small numbers of measurements. 

5.^.1 Velocity derivatives 

Since the ice pack is made up of discrete pieces moving relative to each 
other, the variations in velocity have two types: 

l) for any X, and X, on the same rigid piece, the linear relation 




£S 0 £3 -Stivta \ 




must hold, where oa 


is a constant related to the rate of rotation of the rigid piece. 
2) for and X t on separate pieces. 


u Irt* ) - U-CtC,) 

*+ yv * v 


is 


If we select a random point on an ice floe, 
of velocity. 


arbitrary. 

the partial derivative 


■a;* if+o 8 


is well defined because for small enough % , (?, , and (', almost 

always lie on the same floe and the linearity of the first type of motion 
mentioned above implies the existence of the limit. Thus the partial deriva- 
tive exists, but it describes only the rigid body rotation of a floe. It 
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says nothing about the relative motion between floes which is after all the 
motion causing the opening and closing and the exchange of momentum between 
floes . 

To include the effect of velocity differences between floes we cannot 
let fi become arbitrarily small. Consider the velocity increment 

“ kCk^.-O where is fixed and R begins at 100 km 

and decreases to the smallest value such that (pH-k e , 3) and (x , %) lie 
on separate floes. This increment captures the variatitin in velocity from 
floe to floe and has nothing to say about the rigid rotation of the floe 
(X^y) lies on. The velocity increment will generally decrease as U 
approaches , and we can examine the characteristic rate of decrease by 
evaluating the variance ST ( u(X+£i., - u.(x» }>'&') . If this quantity 

is proportional to for small £ , the velocity at least has a derivative 

in the mean square sense. 

The best data available for examining the variance of velocity increments 
for a range of spacings were extracted from SEASAT synthetic aperture imaging 
radar by Hall and Rothrock, 1981. By comparing two images separated by 
three days, they were able to measure the displacement of enough ice features 
to resolve considerable spatial detail in one dimension of the iee motion. 
After interpolating to evenly spaced points their data have the form 

, VltfJ ] j A* = Z In*- 

Bor 1, the i? -increments were defined as 

Zjt.O- -u.ixd 

v(*ut) - vl *‘) • 

The variance of the increments vas estimated from the 

with the results plotted in figure 1^ . The linearity of the log-log plot 
for small £ supports a power law relationship with vo.r ( X ( 2.) ) r c Pi 1 * with 
-HU Unci }•} is /,3 .It appears 

then that the increments do not decrease as fast as , and therefore 

the velocity does not have a derivative in the mean square sense. 

The statement just made for the variance of increments can be recast in 
terms of the autocorrelation function. It implies that for small r 

6„(r) , /3jW - /-if* 
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thus helping to define the shape of and @_l at r s t> m it r?c , 
and fij_ have the value unity, are continuous, have a continuous first 
derivative of zero, but fail to have a second derivative. 


5.L.2 Velocity differences 

The definition of the derivative as a limit could be discarded in favor 
of a velocity difference over a finite distance, 

u.*x - -L - MtOjA)] 


and similar definitions for V £/L * 

There is no question about the existence of these differences. Defi- 
nitions for strain rate invariants follow naturally by associating U-^v. 
with r ^U-/e>X etc. in the earlier definitions. When L is chosen to be of 
order 100 km, these definitions can be used to describe the large scale 
deformation. 

There are several problems with this approach. First, the L -differences 
carry no information about smaller scales. Second, no basis has been given 
for choosing a particular value for L, . Third, measurements of velocity are 
rarely available at uniformly spaced points. Still most analyses of sea ice 
strain have been done with a length scale implicitly fixed by the scale of 
the observations. To interpret these analyses we must determine the role 
played by L ■ 

Consider the covariances between the various L -differences. These 
follow from the covariance structure of the velocity itself. For example 






+ 0)1 ( r - L )j . (5) 


To interpret this equation, suppose the J_ —difference is measured 

at two points separated by a distance T along the pt axis. The right hand 
side of equation S approximates the second derivative of Gjj at f . 

For small T , t3| is concave down so the two L -differences are positively 

correlated. At some value of T , Qq (r) has an inflection point. Over 
such a distance the L -differences become uncorrelated; at longer distances 
the -differences are negatively correlated. This example has been worked 
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for the V>x,L component at points separated along the pL -axis, but the 
same procedure can be used to correlate any two first velocity differences 
at any two points. 

It is apparent from this and figure M> that quantities related 

to the /. -differences are highly correlated in space only over distances 
of less than about LOO km. Furthermore this statement is not sensitive to 
the choice of £ , since the expression S will change sign at approximately 
the same ■T , independent of L . The correlation length scale for velocity 
by comparison is roughly twice as large (see figure l£)» 

The ice velocities are well correlated at distances of up to 1000 "km, - 
as evidenced by the patterns in the long term ice motion (figures 2-h ) and 
in shorter term motions (see Thorndike and Colony, 1980 for instance). It 
may seem surpising not to find similar patterns in the strain fields. The 
reason why such patterns are not found Is that the strain field has, by 
equation 5, and figure 16 , a correlation length scale of only a few hundred 
kilometers. Measurement arrays like the 1979 and 1980 buoy arrays with 
a buoy spacing of roughly 500 kilometers are too coarse to resolve patterns 
in the strain field. 

For the example just worked we can find the variance of the L -difference 
by setting - 



The symmetric matrix of covariances between pairs of L -differences at the same 
point is 






•f [l-6„(L)J 


Vfri. 




^[l +Wfa0-2B(0] 


O 




l[) f 6(<3L)Ve(lv>] 


V 

£3 = ( % +" Q-) • 

This matrix can be evaluated at L — L00 
from the plotted functions (?jj and @J_ 


$ [i- a.wj 

km for example by reading points 
in figure 1^. 


2 ? 


ORIGINAL PAGE IB 
OF POOR QUALITY 



Moo W 

- a i -2. 
pi =• //> «& 


The covariance matrix for the L -differences illustrates several properties 
of the way the ice deforms. Notice that the variances of the parallel dif- 
ferences U. Xj L and are smaller than the perpendicular differences and 

Vjc # j_ . The parallel differences relate to the linear stretching and 
contracting of the ice pack. The perpendicular differences relate to shear 
and rotation. Apparently the second type of motion is greater, "but recall 
the word of caution given earlier that the interpretation of this kinematic 
evidence is ambiguous. We cannot say whether the observed motion reflects 
some property of the ice pack or some characteristic of the driving forces. 

With the L -differences one could define [_ —strains in a natural way. 

The statistics in Table XV of the -strains follow from the above matrix. 

These values predict, on the basis of the spatial correlations above, that 
the large scale divergence will typically be smaller than the vorticity and 
shear, a prediction verified by many different sets of observations. 


TABLE IV 

Topical strain rates as a function of gauge length L. . 

L-sVi.V, L-sW in ftrceJ ycr V*rK*a 
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5.4.3 Deformation measurements 

Several studies have "been made over the past decade or so using obser- 
vations of velocity at an array of stations, with the objective of monitoring the 
large scale deformation of the region covered by the array. While it is never 
stated just this way, the implicit definition of large scale deformation appears 
to involve spatially averaged derivatives over the region. 



R 


where A is the area of the region R sampled by the array. Similar expressions 
define the other large scale average velocity gradients. Application of the 
Green-GausS theorem implies 

& * i: / “• m-l & ^ 

C 

where *W is the outward normal to the perimeter C of ft , and % a unit vector 
in the ^ direction. 

Note the difference between the L -difference and the large scale 
average derivative. The former can be measured exactly as ufay-i, f ) - ; 

the latter requires measurements at every point around a closed curve. In the 
next section we will discuss how many measurements are required to achieve a 
desired accuracy in the large scale deformations. But first we review some of 
the results of the studies just alluded to. 

In these studies the deformation estimates were made by finding the 
linearly varying velocity field which most closely fit the observed velocities 
at an array of points within R, 

u(*)= u 6 + 

where the matrix M contains the four large scale average velocity partial deriva- 
tives. M and U 0 are determined using a least squares fit of the observed 
velocities u and positions % . The strain rate invariants were then computed 
from the elements in ff . 

The values, given in Table V, confirm some of the results deduced above 
from the observed spatial correlation functions. For example the standard 
deviations of vorticity and shear are generally larger than the standard 
deviation of divergence. Note also that the mean divergence is small com- 
pared to the mean vorticity and shear. 
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TABLE V 


Statistics of Sea Ice Deformation. 

All deformations are expressed in per cent per day, (IjS/d — I.l6xi0” 7 s” 1 ) 


Length scale 

diameter 

(km) 

Season 

Divergence 

Shear 


Vorticity. 

* .. 

Strain 

magnitude 

mean 

r- 

mean 

st.dev. 

S' * 

mean 

st.dev.’ 

mean 

st.dev. 

.20 

spring 

— 

0.65 

— 

0.86 

— 

1.1 

— 

100 

spring 

— 

0.43 

— 

0.54 

— 

0.61 

— 

200 

| winter 

0.07 

1.0 

1.6 

1.6 

-0.52 

2.0 

2.1 


[ summer 

-0.03 

1.6 

3.5 

2.2 

— 

— 

— 

800 j 

f winter 

0.02 

0.56 

1.0 

0.86 

— 

— 

1.3 

i 

1 

L summer 

0.16 

1.0 

1.6 

0.86 

— 

— 

2.2 

800 ) 

' winter 

0.07 

0.58 

0.77 

0.54 

— 

— 

1.1 

i 

summer 

0.06 

0.80 

1.0 

0.68 

— 

— 

1.5 

1000 

winter 

0.02 

0.45 

0.95 

0.82 

-0.78 

1.5 

1*1 


summer 

-0.11 

0.76 

1.5 

0.95 



1.9 


Experiment 


AIDJEX 1972 

AIDJEX 1972 

1975 AIDJEX 
manned array 


1975 AIDJEX 
Buoy array 


1962-64 T3, 
irp-10, 
Arlis II 


1979 buoys 


Sources: Hibler, et al. , 1974; Colony and Thorndike, 1981; Maykut, 1981. 
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The distribution of the invariant & is plotted in figure 20 . Note 
that nearly always % < 8 < . In this range , the two principal values 

of the strain rate tensor have opposite signs. The figure indicates that 
it is almost always the case that when the ice pack is extending in one 
direction, it is contracting in the perpendicular direction. Hibler et al. , 
197 ^ also comments on this for strain measurements on a 20 kilometer scale. 

It is also evident from Table IT that the summer values are usually 
somewhat greater than the winter-spring values. Presumably the ice pack 
is weaker and offers less resistance to deformation in the summer. (The 
alternative hypothesis, that the driving forces are larger in the summer, 
is not true. ) 

The mean quantities in Table -V produce large strains over the course 
of a year. The year long deformations at the AIDJEX arrays are shown in 
figure 2./ . The region experienced a net clockwise rotation of about 35° • 

The principal deformation involved a stretching of about 90$ in the east-west 
direction and a contraction of about ^- 0 % in the north-south direction. The 
net divergence was not significantly different from zero. The two nested 
arrays experienced similar deformations. 

Deformation estimates at a number of points in the central basin for 
the years 1979 and. I 98 O are presented in figure 22. The strain ellipses 
typically show large shear and small divergence. The only pattern evident 
in the figure is the similar alignment of the major principal axis of shear 
for the five points closest to the pole. Perhaps we should not expect to 
see any patterns in the deformation displayed on this scale. Recall that 
the spatial correlation function for ]„ -differences of daily velocities 
has a length scale of only about 400 kilometers. The deformations at points 
separated by greater distances should evolve more or less independently. 
Conceivably there may be some underlying spatial pattern to the deformation 
rate which is small on a daily basis but which accumulates to a large spa- 
tially organized strain over the course of a year. The strain ellipses 
for 1979 and 1980 do not support this idea. On the other hand the long 
term average vorticity is clockwise throughout the central hasin. It amounts 
to about L0° per year. 

5-Av^ Interpretation of deformation measurements 


A number of authors have confronted the difficulties oi describing the 
deformation of this dis continuum. Nye ( 1973 ) for instance, puts forward a 
definition of "strain on a length scaled By first smoothing the velocity 
field using a kernel of length scale lu 9 he obtains a new velocity field 
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which is differentiable and for ■which the usual notions of deformation "based 
on partial derivatives are valid. Papers "by Hibler et al, , 197^ and Thorndike 
and Colony, 1977 take a similar point of view "by attempting to partition the 
ice deformation into two parts. The first part is assumed to be associated 
with the spatial variation of the atmospheric and oceanic forcing fields and 
with the geometry of the basin, to have a characteristic length of 10 km, 
and to be differentiable. The second part is associated with the irregular, 
discrete geometry of the ice pack. It has length scales less than 10 km and 
is characterized by discontinuous variations in velocity. In these papers, 
the first kind of variation is referred to as the large scale, underlying, 
continuum, or linear part of the velocity field, while the second is referred 
to as the small scale, local, perturbation, fluctuating, or non-linear part. 

Velocity variations of the first type are regarded as signal and those of the 
second type as noise. This nomenclature reflects the hope that any physical 
process of interest can be parameterized in terms of the large scale signal, 
with small scale noise only making it difficult to measure the signal well. 

This convention obscures the true nature of the velocity field which 
is that it has variations on all length scales, with a smooth decrease in 
amplitude for decreasing length scales. There is no clear division between 
large and small scale. Furthermore there is no clear reason to associate 
large scale with signal and small scale with noise. For studies involving 
the actual opening and closing of leads, the small scale phenomena may indeed 
he the signal. Further we must be prepared for the possibility that the small 
scale phenomena cannot he parameterized in terms of the large scale motion. 

Perhaps a better conceptual model of the spatial structure of the ice 
velocity is a system which accepts a smooth input and produces a discontinuous 
output. simple example of such a system is the greatest integer function: 
jjxJ = greatest integer less than or equal to 
The shift in emphasis from the earlier model is this. The earlier model 
viewed the velocity as the sum of separate contributions, one smooth, one 
discontinuous. The alternative views the velocity as a discontinuous response 
to a smooth input. This point of view may lead more naturally to phenome- 
nological descriptions of the properties of the ice pack, through a comparison 
of the input and output fields. To my knowledge this has not been attempted 
because data with adequate spatial resolution of the velocity field are still 
too scarce. Such a study would not explain why the velocity field had 
certain properties. That explanation must be based on rather deeper understanding 
of the geometry of the ice pack and the forces which act between floes than we 
have at present. Still it would he useful to compare some of the properties 
of the ice velocity with properties of the external driving forces, i.e., the 
geostrophic wind and the ocean currents. 
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5.4.5“ . Errors in estimating the large scale deformation 

The large scale deformation estimates of the previous section are 
subject to sampling variations since they depend on the particular array 
of measurements points. There have been several attempts to quantify the 
sampling errors. Hibler et al. , 1974; Wye, 1975; and Thorndike and Colony, 

1977 examined the departures of the measured velocities from the best fit 
linear relationships. These departures, called residuals, homogeneity 
variations , or nonlinear fluctuations , . represent the variability of the 
actual velocity field on scales smaller than the diameter of the region of 
interest. 

Prom observations over a 20 km scale taken during the spring of 1972 > 

“1 

Hibler et al. found root-mean- square velocity residuals of 0.06 cm sec , 

compared to typical linear changes in velocity ( ) over 20 km of 

“1 

O.lk cm sec * Thorndike and Colony used observations from a 100 km scale 

—1 

taken during the spring of 1975 and found rms residuals of 0.4 cm sec and 

typical linear changes over 100 km of 1*1 cm sec \ Their summer values -were 

—1 —1 

somewhat larger: 1.1 cm sec and 1.8 cm sec for the nonlinear and linear 

contributions. These residuals can he regarded as errors with respect to the 
large scale average derivatives for the region. When the number of observa- 
tions is small, the estimated large scale deformations are strongly contaminated 
by these errors. 

Using what we know about the correlation functions for velocity, we can 
profitably address the sampling question from a different angle. Taking the 
line integral definition for the large scale velocity derivatives , we ask: 
how many points around the perimeter of the region must he sampled to resolve 
the integral to some desired accuracy? Intuition suggests that the measure- 
ments should be spaced closely enough to permit good interpolation hut not 
so closely as to he highly redundant. A correlation between velocities at 
neighboring measurement points of 0.5 might be a fair guess; this would 
indicate a spacing of about 400 km. 

To get a better answer, we estimate the line integral in equation 6 
around a circular region of radius r as 


0), 


M 


M 


U-l 


B: 


fc- 


L * • V M • 




Then as the number of measurements H increases, approaches 

Since "3^ is a linear combination of velocity components, we can find its 
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variance in terms of the variances and covariances of the velocity components 
at all M points, and these can he evaluated in terms of the correlation 
functions 8„ and . Let U. he the vector involving all 2M velocity 
components, A he the vector of coefficients, and ft. he EILUj . Then 
D h =A T y. maEoJv P? ft. A • This calculation has heen done for 


M" io which was large enough to resolve ~b u. /~& 7 C . The dependence of 
on the radius *r is shown in figure 23a. 

This curve depends on 6 ( j and which for this purpose were taken to he 
(l-r/ft.) with h- = 2000 km, ia = 1300 km for and d = 900 km 

k = 1000 km for fij, * These functions incorporate the features of the 
observed correlations in figure 16. The leading factor ( ] — r/a. ) makes 
the assumed functions vary linearly near r = 6 . The behavior of the 

correlation functions near is related to the variance of increments 




Vo.r ( u(r+Zi) - u(t^) * 5L^O- e i|(4>) 


* 2^ A/o, for sma11 A 

■which is a rough approximation to figure 19 "where the variance of increments 
is seen to he proportional to . The results in figure 23 are quite 
sensitive to the behavior of the correlation functions. They should not 
he taken too seriously until more is known about Gj| and G[_ near A~Q . 

With some caution. then we interpret figure 23a as follows. Typical 
values for the large scale average derivatives are X% per day decreasing 
somewhat with the size of the region over which the average is taken. The 
derivatives and are generally larger than 

"fcu./-** and ~iv/ as a consequence of 6j_ falling off more 

rapidly than fijj . 

The variance of the error in estimating using only M points, 

E ( D M - t^o) can he calculated in the same manner simply hy redefining 
the coefficients in the vector A . These results are presented as a 
fraction of the signal variance E(D^) in figure 23b. When M is 3, 
as was the case for AIDJEX 1972, the 1975-76 manned AIDJEX array, and 
LOBEX 79 , the ratio of estimation error variance to signal variance was 
about 0.7. 

With six measurements around the circumference of the region, the 
average derivations are resolved fairly well, and only small improvement 
is realized by adding more measurements. The results are insensitive to 



Suppose a region of interest is intersected by a number of cracks, each 
of 'which is opening or closing at some rate. If the rates were known, the 
total rates of opening and closing for the region could be found. In practice 
we are not able to measure the motion at each crack, but only the motion of a 
few points in the region. How are we to use these few measurements to esti- 
mate the opening and closing? 

In an analogous situation in one dimension, we have k. cracks each with 
opening or closing rate UL^ L= • * * , k, . We can imagine the velocity U. to 
be a random function of pi having random discontinuities at random points 
and being constant between the points. Suppose we have measured only the 
motion at the end points of the region, an interval of length L . Then 






(“•L , o) 

(w-i, O 




A + t-- X u i , 


and the problem is to estimate A and C given U . At first sight the 
situation seems hopeless. LL clearly contains information about the net open- 
ing or closing but not about the total opening and closing. However, knowledge 
about the random variables )? U.,,^*** can be used to. make probabalistie 
statements about the opening and closing. Suppose for instance k has the 
Poisson distribution with parameter A . This means that is the average 

number of cracks in an interval of length L . Then the probability of 
finding ^ cracks in a random interval of length )_ is 
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If AL = 0 . 1 . for instance, the probability of getting ft= 0 is 0.9051, 

= 0.0905, and prt>L(,k>i) = O.O 0 L 5 . Therefore, with high probability, 
there is either no crack or just one in the interval and the observed value 
of the velocity difference UL -would itself indicate the total opening and 
total closing. Of course for larger it becomes more likely that several 
cracks intersect the test interval, in which case the observed U. cannot 
separate the opening from the closing. 

A similar approach for the two dimensional problem is to imagine the ice 
pack to be crisscrossed by a family of random cracks, defined by the random 
straight lines 

Pi + jj Jwfy “ 

where each has the uniform distribution on (e>, JLy) and the normal distances 

from the origin to each line form a sequence of Poisson points. Imagine 
that associated with each crack is a velocity discontinuity U- b having 
the Gaussian distribution. This describes a random vector field having 
structure akin to the velocity field of sea ice. The observations of Hall 
and Rothrock can be used to estimate the parameters in the Poisson and 
Gaussian distributions. These random fields have the following properties. 

They consist o,f discrete rigid floes. The floes have a distribution of sizes 
determined by the Poisson field of lines. The Poisson lines are isotropic 
and homogeneous. The velocity difference between any two points is the 
vector sum of the velocity discontinuities encountered getting from one point 
to the other. 

At each crack the opening or closing is determined by the projection 
of the velocity difference vector onto the normal to the crack. 

opening = max ^ d, u.^ ) 

closing = min ( o, u;Ci»0l al ~ «. &c ) 

Thus it is a simple matter to evaluate the total opening and closing for any 
realization of the random field. 

Of course one can also imagine measuring the velocity at a few points 
and constructing the -strain rate invariants from the observed velocities. 

In this way one can test for a relationship between the -strain rate 
invariants calculated from a few sampled velocities and the total opening 
and total closing found by tallying up the activity at every crack. 
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In an attempt to carry out this program, parameter values were taken 
to he /\ =D.0?km -1 , ^ = 0, <T 2 = 1 (arbitrary units) 2 . An L = 100 km 
triangle was used to simulate the sampling procedure used during AIDJEX. 

A large number of random fields were generated. 

Each realization of the random field is defined by the sequences of 
random numbers k. } &i , U. • , # Tj; for i " ij ■ • > k . Here the &i are drawn 

independently from the uniform distribution on ( b >ir) . The values for 
i'l form a Poisson process with parameter /l . This is achieved by drawing 
the increments 4- -T* t independently from the exponential distribution 

with density . The process is terminated as soon as exceeds 

100 kilometers since none of the subsequent .lines would intersect the 100 
kilometer region. Finally the and were drawn independently from 
the normal distribution with zero mean and unit variance. 

For each field the velocity was measured at the three specified points 
and the L,-strain rate invariant s |£| and 6 were evaluated. Also for 
each field the total opening and closing were evaluated using eq. T» 

These were normalized by |£1 and displayed versus 9 in figure 25 . For each 
realization of this random field two points are plotted, ( B , A/m) and 

G/jEl ) . From the plot it is clear that there are not unique values 

of A/iE) and C/lE'l corresponding to a given 0 . Instead there is a 
distribution of values for A/]E\ (and f or C /\E) ), and this distribution 
changes with 0 . The distribution^ sketched in figure 2$ ; are broad in the 
sense that probable departures from the mean are at least as large as the mean 
The interpretation of this exercise is that the total opening and closing 
are only weakly determined by the L -strain rate invariants based on three 
measurements 100 kilometers apart. Had the relationship been a strong one 
the distributions in figure 25 would have been narrower, or to say it differ- 
ently, the points in figure 25 would have clustered more closely around 

curves like those in figure /? . It may still be useful to imagine these 

smooth curves but only with the recognition that the actual opening and 

closing scatter widely about the imagined curve: 


fe) + random error 

C/ , , fe \ + random error 

/IF} ' 


where the random terms and the B - dependent terms make roughly equal contri- 
butions to the total opening and closing. 
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Some of the scatter in these figures arises from having used only 
three points to estimate the deformation invariants. With more than three 
points, the errors in the invariants could "be reduced, somewhat as shown 
in figure 23. But even if the invariants are determined without error, 
scatter will remain in figures 2k and 25. This is because many different 
fields of motion with different total opening and closing could have exactly 
the same large scale average deformation. 

The total opening and closing are an essential part of the theory of 
the ice thickness distribution. The above results suggest that it may not 
be satisfactory to represent the total opening and closing as functions 
of the large scale average velocity derivatives. When many velocity measure- 
ments are available in a region, it is possible to estimate more of the 
spatial structure than simply the average velocity derivative. An attempt 
to estimate the variability of velocity within the region may help to 
reduce the scatter which is anticipated on the basis of figure 2k, Or 
perhaps it will be enough to assume that the variability is constant in 
time. One might then drive the thickness distribution calculations with 
opening and closing time series which have the right statistical properties 
even if they may have large errors on a day to day basis. The statistical 
properties can be inferred from figure 2k, Of course this figure is the 
consequence of a particular conceptual model motivated by a limited data 
set — the SEASAT SAR data from early October 1978* More data need to be 
collected and studied before these ideas can be extended to other times and 
places with confidence. 

5.^.^ The relationship between kinematics and stress 

Kinematic data can be used for studying the relationship between 
stress and strain for sea ice. Suppose for instance that the state of stress 
tensor CT suitably averaged over some region, is related to kinematic quan- 
tities e and the ice geometry S by some expression of the form 

(T =• F(£,s) . 

Although <r cannot be measured directly, its divergence V-C can be inferred 
indirectly from the observed momentum balance when all the other terms in 
that balance are known. If a particular function F is hypothesized, it 
can be tested by evaluating both sides of the equation ^ f r ( £ J s ) , 

the left hand side as a residual from the momentum balance and the right 
hand side in terms of observed kinematic and ice geometry quantities. This 
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was one of the objectives of AIDJEX (Maykut et al. , 1972). In practice 
the test is difficult to make since neither V-(T nor V-F can he deter- 

mened very accurately from observations. See Bothrock et al., 1980 . 

Under special circumstances, meaningful tests may be possible, however. 
Uhen the ice deformation is strongly divergent, ice floes tend to move apart. 
With no floe-to-floe contacts there can be no S7- <T forces. The left hand 
side should differ from zero only by measurement error. These errors are 
probably small enough to permit a useful test. 

Another special situation of interest is when the ice is being forced 
up against the coasts by the wind. As the ice moves toward the coast it 
converges, becomes stronger, and eventually becomes strong enough to resist 
further deformation. If the on shore winds persist, a zone of motionless 
ice can widen to several hundred kilometers (Pritchard, 1977)* In this situa- 
tion, intuition suggests that the V*(T vector should point off shore and 
should increase in time until it balances the on shore wind stress. With 
simultaneous ice trajectories at 100 , 200 , 300 and 400 km from the coast, 
it should be possible to observe the amount of convergence required to pro- 
duce the required resistance to the wind. 

At greater distances from the coasts, there is evidence that the ice 
stresses embodied in V*<T are usually small. Little success is anticipated 
in trying to observe them. Any stress-strain law which provides adequate 
resistance to deformation near the shore is probably adequate for full basin 
dynamic modeling. The observations may never be adequate to discriminate 
very selectively between candidate stress-strain laws. 

3.5 Discussion 

I have tried in this paper to bring together some of what is known 
about the motion of sea ice. The emphasis has been on the departures of 
the true motion from the long term mean circulation. This is a compari- 
tively new topic, the investigation of which has only recently been made 
possible, mainly by the developments of satellite positioning techniques. 

Perhaps because of its newness, the field lacks a clear agenda of questions 
to be answered. Instead we are still trying to characterize the motion, to 
determine the magnitude of the motion on different scales and to identify 
motions which have some intrinsic interest or are related to other phenomena 
of interest. 

The time and space correlations have been used extensively here for 
several reasons. First they give a compact description of the motion. In 
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principle, the correlation function involves one time and two space variables, 
"but in practice a great deal of information can he extracted from three 
functions, each of a single variable: 0^[ t ) and . Second, these 

three functions are accessible to observation. Third, properties of many 
kinematic quantities can be deduced from these functions, as illustrated 
in the text. Fourth, these functions form the rigorous basis for answering 
questions related to interpolation, prediction, and experiment design. 

Many of the results presented here are based on sample autocorrelation 
functions deduced from limited data. The data available for estimating 
correlations at small space lags (less than 100 km) are meagre indeed. 

More work along the lines of Hall and Rothrock would help to resolve this 
part of the correlation function. The behavior of the correlation functions 
0l|(r) ckwi (?j L (-r) in the limit of small r is an important property of 
the motion related to the granular nature of the ice pack. 

As mentioned in the text, care must be exercised in choosing correlation 
functions or contradictions (negative variances) can occur. In fact this 
has happened in some of the calculations done using the correlations tabulated 
in Table with linear interpolation to intermediate distances. This means 
that the piecewise linear function defined in the table is not positive 
definite. A useful objective would be to find a positive definite analytical 
form which approximates the observed correlations including the behavior 
as 1* approaches zero. 

The Pois son-Gauss model, presented as a way to study the relationships 
between the local opening and closing and the large scale deformation, has 
suggested that the relationship may be weak. Observations against which 
to test this suggestion are sorely needed. 


Table VI. Spatial correlation functions for sea ice velocity. 


Distance (km) 


6a 




0 

1.00 

1.00 

100 

• 98 

.95 

200 

.91 

.81* 

1*00 

.68 

.51 

800 

.37 

.06 

1200 

.19 

-.09 

1600 

.10 

-.10 

2000 

.01 

-.06 

2 l *00 

.00 

.00 


38 


Acknowledgment . This work was supported by the National Oceanic and 
Atmospheric Administration Grant NA80-AA-D-00015 , which was funded in part 
"by the Global Atmospheric Research Program and the Office of Climate Dynamics, 
Division of Atmospheric Sciences and the Meteorology Program, Division of 
Polar Programs, of the National Science Foundation, and the Office of 
Naval Research, Arctic Programs. The paper also reports on work supported 
by the Office of Naval Research Contract N0001-4-79-C-0^l8 . 

Some of the ideas discussed here grew out of informal discussions with 
Drew Rothrock and with Roger Colony. I take this opportunity to acknowledge 
their contribution? and to thank them. 




ORIGINAL PAGE SS 
OF POOR QUALITY 

References 

AIDJEX Staff, 1972. Station positions, azimuths, weather, 1972 AIDJEX Pilot 
Study, Preliminary data, AIDJEX Bulletin go. lU, 63-208. 

Colony, R. and A. S. Thorndike, 1980. The horizontal coherency of the motion 
of summer arctic sea ice. J. of'Phys, Oceanography , 10(8), 1281-1289* 

Colony, R. and A. S. Thorndike, 1981. Sea ice strains during 1979* Proceedings 
of the 6th International Conference on Port and Ocean Engineering Under 
Arctic Conditions, Quebec City, Quebec, 27-31 July 1981, 619-628. 

Craig, L. W. , 1972. High frequency accelerations of arctic pack ice. Master's 
• thesis. Department of Civil Engineering, University of Washington, Seattle. 

3 

I Doronin, Yu. P. and D. E. Kheisin, 1975* Morskoi Led. (English translation, 

i Sea ice, 1977* Amerind Publishing Co. Pvt. Ltd., New Delhi.) 323 pp. 

Hall, R. T., 1980. AIDJEX modeling group studies involving remote sensing 
data. Sea Ice Processes and Models , R. S. Pritchard, Editor; University 
j of Washington Press, Seattle. 

j Hall, R. T. and D. A. Rothrock, 1981. Sea ice displacement from SEASAT 

synthetic aperture radar. J. of Geophys, Res. . 86(CH) , 11078-82. 

! 1 

j Hastings, A. D. Jr., 1971. Surface climate of the Arctic Basin. Report ETL- j j 

TR-71-5, Earth Sciences Division, Geographic Sciences Laboratory, U.S. Army 
Engineer Topographic Labs, Fort Bel voir, VA. 

Hibler, W. D. Ill, W. F. Weeks, A. Kovacs, and S, F. Ackley, 1974. Differential 
sea-ice drift. I. Spatial and temporal variations in sea ice deformation. 

J. of Glaciology , 13(69), 437-455. 

Hunkins, K. , 1967. Inertial oscillations of Fletcher's Ice Island (T-3). 

J. of Geophys. Res. . J2(4), 1165-73. 

Hunkins, K. and J. Hall, 1971. The navigation of Fletcher's Ice Island (T-3). 

AIDJEX Bulletin No. 7 . 47-58. 

Hunkins, K. T. Manley, and W. Tiemann, 197 9* Observations of position, ocean 
depth, ice rotation, magnetic declination and gravity taken at the Fram I 
drifting ice station. CU-1-79 Technical Report No. 1 . Office of Naval 
Research, Department of the Navy, 60 pp. 

Kolmogoroff, A. s 1941. The local structure of turbulence in incompressible 
viscous fluid for very large Reynolds numbers. Comptes Rendus (Doklady) 






References, continued 


ORIGINAL PAGS 19 
OF POOR QUALITY 


Martin, P., C. R. Gillespie, A. S. Thorndike, D. Wells , 1978. Position 

measurements of AIDJEX manned camps using the Navy Navigation Satellite 
System. AIDJEX Bulletin No. 40 , 83-102. 

Maykut, G. A., 19 81. Large scale heat exchange and ice production in the 
central arctic. To appear in J. of Geophys, Res. 

Maykut, G. A., A. S. Thorndike, and N. Untersteiner, 1972. AIDJEX scientific 
plan. AIDJEX Bulletin No. 15 , 67 pp. 

McPhee, M. G. , 1978. A simulation of inertial oscillations in drifting pack 
ice. Dynamics of Atmospheres and Oceans , 107-122. 

Nye, J. F. , 1973* The meaning of two-dimensional strain-rate in a floating 
ice cover. AIDJEX Bulletin No. 21, 9-17. 

Nye, J. F., 1975. The use of ERTS photographs to measure the movement and 
deformation of sea ice. J, of Glaciology , 73(15? , 429-436. 

Nye, J. F. , 1975. A test of the ice thickness distribution equations by 
measurements of ERTS pictures. AIDJEX Bulletin No. 28 , l4l— 150. 

Nye, J. F.-, 1976. A coordinate system for two dimensional stress and strain 
rate and its application to the deformation of sea ice. AIDJEX Bulletin 
No. 33 , 131-143. 

Popelar, J., J. Kouba, and D. Wells, 1981. LOKEX Satellite Positioning, 

Draft report. Earth Physics Branch, Department of Energy, Mines and 
Resources, Ottawa, Ontario, Canada. 

Pritchard, R. S., 1977. An estimate of the strength of arctic pack ice. 

AIDJEX Bulletin No. 34 . 94-113- 

Rothrock, D. A., 1979. Modeling sea ice features and processes. J. of 
Glaciology , 24(90), 359-375. 

Rothrock, D. A., R. Colony, and A. S. Thorndike, 1980. Testing pack ice 

constitutive laws with stress divergence measurements. Sea Ice Processes 
and Models , R. S. Pritchard, Editor', University of Washington Press, 102-112. 

Rothrock, D. A. and R. T. Hall, 1975. Testing the redistribution of sea ice 
thickness from ERTS photographs. AIDJEX Bulletin No. 29. 1-20. 

Sverdrup, H. U. , 1926. Dynamics of tides on the north Siberian shelf, results . 
from the Maud expedition. Geofysiske Publikajoner IV(5)» 75 PP- 

Tabata, T. , T. Kawamura, and M. Aota, 1980. Divergence and rotation of an 

ice field off Okhotsk Sea Coast of Hokkaido. Sea Ice Processes and Models, 


R. S. Pritchard, Editor; University of Washington Press, 273-282. 


References , continued 


ORJGfNAL PAGE IS 
OF POOR QUALITY 


Thorndike, A. S., 1973. An integrated system for measuring sea ice motions. 

In Ocean 73 . IEEE International Conference on Engineering in the Ocean 
Environment, IEEE Publ. No. 73 CHO 774-0, 490-499. 

Thorndike, A. S. and J. Y. Cheung, 1977* AIDJEX measurements of sea ice 
motion 11 April 1975 to May 1976. AIDJEX Bulletin No. 35 , 1-149. 

Thorndike, A. S. and J. Y. Cheung, 1977- Measurements of sea ice motion 
14 May 1976 to 30 November 1976, October 1975 to December 1976, and 
January 1977 to September 1977. Unpublished data reports, available 
from the authors. 

Thorndike, A. ,S. and R. Colony, 1977. Estimating the deformation of sea ice. 

In Proceedings of the Fourth Annual Conference on Port and Ocean Engineering 
under Arctic Conditions , Vol. 1, D. B. Muggeridge, ed. , 506-517* Memorial 
University of Newfoundland, St. John’s. 

Thorndike, A. S. and R. Colony, 1980. Large scale ice motion in the Beaufort 
Sea during AIDJEX, April 1975 to April 1976. In Sea Ice Processes and 
Models , R .S . Pritchard, ed.. University of Washington Press, 249-260, 
Thorndike, A. S. and R. Colony, 1980. Arctic Ocean Buoy Program, data report, 

19 January 1979-31 December 1979. University of Washington, 131 pp. 
Thorndike, A. S. and R. Colony, I98I. Arctic Ocean Buoy Program, data report. 

1 January 1980-31 December 1980. University of Washington, 127 PP* 
Thorndike, A. S. and R. Colony, 1982. Sea ice motion in response to geostrophic 
winds. To appear in J. of Geophys. Res, 

Thorndike, A. S., D. A. Rothrock, G. A. Maykut and R. Colony, 1975* The 

thickness distribution of sea ice. J. of Geophys. Res. , 80,(33), 4501-4513. 
Tucker, W. B. , W. F. Weeks, A. Kovacs, A. J. Gow, 1980. Nearshore ice motion 
at Prudhoe Bay, Alaska. In Sea Ice Processes and Models . R. S. Pritchard, 
ed.. University of Washington Press, 261-272. 

Zubov, N. N. , 1943. L'dy arktiki (Arctic ice), Glavsev-morput (Northern Sea 
Route Administration), Moscow. Translation 1963. U.S. Navy Hydrographic 
Office Translation 217, available as AD-426-972, from NTIS, Springfield, VA. 



p, 1 

Figure 2. 
» 

Figure 3. 


Figure 
Figure 5- 


Figure 6. 


Figure 7. 
Figure 8. 

Figure 9- 


Figure 10. . 


Figure 11. 
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The general circulation of ice in the Arctic Ocean (from Doronin 
and Kheisin). 

Buoy trajectories in 1979. She numbers indicate the months in which 
each buoy was deployed and failed. Trajectories ending in an arrow 
continued to operate in 1980. Plotted points correspond to the beginning 
of each month. 

Trajectories for drifting buoys during a) 1980, and b) 1981. 

(a) Latitude of ice station Caribou, versus time, 2 Sept-l6 Hov 1975. 

HavSat data, sampling interval, 2.5 days. 

(b) Latitude of ice station Big Bear, versus time,- 2-5 Sept’ 1975* 

HavSat data, sampling interval, U hours. 

(c) Arbitrary y-coordinate of Big Bear, versus time, 2-3 Sept 1975. 

Acoustic tracking data, sampling interval 15 minutes. 

(d) Arbitrary x-coordinate of Big Bear, versus time, 2 Sept 1975. 

Acoustic tracking data, sampling interval one minute. 

Only in (d) is the measurement error visible (±3 m) . 

(a) The variance of the increments and - Jftb) for 

ice station Caribou 1975-76. yt and h are arbitrary Cartesian 
coordinates. ® 

(b) As in (a) using precise acoustic tracking at ice station Big Bear, 
an 8-day period, late summer 1976. 

Velocity histograms (from Thorndike and Colony, 1980). 

The power spectral density of the velocity of ice station Caribou 1975-76. 

The units of spectral density are velocity a /frequency, cm* sec -1 . The 
total velocity variance for these data is 1^5 cm* see*" 3 . 

The power spectral densities of the geostrophie wind and the ice velocity, 
from drifting buoy 1901 during 1979. The units of spectral density are 
cm sec-' . The variance of the wind is 53 m 2 see -2 . The variance of 
the ice is 59 cm 2 sec -2 . 

The power spectral density of the velocity of ice station Big Bear 

during summer 1975. The variance is 173 cm 2 sec -2 . Units of spectral 

density are m 8 * s~ f . (Linear frequency scale, from Colony and Thorndike, 1980). 

The complex time correlation function for ice station Caribou, 1975-76. 

The velocity variance is 1^5 cm 2 sec” 2 . 


Figure 12. The time correlation function for the u. and V components of velocity 

from 28 grid points in the central part of the Arctic Basin, from drifting 
buoy data collected during 1979. The cross correlation between u. and 
V is negligible. Variance of M.= 23 cm* sec -5 ; variance of v = 22 cm 2 sec -2 


Figure 13. (a) The real and imaginary parts of the weight to be used for inter- 

polating to the mid point of a time interval of duration V . The ratio 
of the measurement error to the standard deviation of velocity is <r/p. . 

(b) The variance of the interpolation error expressed as a fraction of 
the variance of the ice velocity. Data from figure 11 were used here. 
Logarithmic time scale. 


ORIGINAL PAGE ffi 
OF POOR QUALITY 


Figure lU. 
Figure 15. 

Figure 16 . 
Figure 17. 

Figure 18. 

Figure 19. 
Figure 20. 

Figure 21. 

Figure 22. 
Figure 23. 

Figure 2k. 
Figure 25. 
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(a) The weight for prediction over time -r . 

(b) The prediction, error variance versus prediction time. See Figure 13 
also. 

A schematic representation of the ?C and $ components of velocity at 
the points (o,e> ), (f,D), (e» ,■*-). Isotropy implies £ac = EL-?, 

£ U j " Eae. , and "ELe. = E Jf , but does not imply equality between 

E-fcc andEae or & If- and E bi 

Observed correlations between velocity components as functions of distance. 

The variance of the interpolation error as a function of the grid spacing, 

J_, and the measurement error, c r* . 

The strain rate invariant 0 indicates whether the motion is dominated 
by divergence, shear, or convergence. The total opening and closing 
are sometimes treated as functions of 0 . 

The variance of velocity increments versus interval length, from SeaSat 
synthetic aperture radar data, 3-5 October 1978, Beaufort Sea. 

Histograms of the strain rate invariant & . Solid line is for drifting 
buoy data from 1979, grid spacing about 500 tan (Colony and Thorndike, 1981). 
Dots are for AIDJEX manned camp data, spacing 100 km (Nye, 1976). v/jift: 2.6 

Strain ellipses for the AIDJEX buoy array, 800 tan diameter (a); and for 
the AIDJEX manned camp array, 200 tan diameter (b). The ellipses show 
the deformation of a circular region on 1 May 1975, to the date indicated. 

The principal axes of the ellipse are the principal strains. The angle 
from the horizontal broken line to the major axis is the principal 
direction. The rigid body rotation is indicated by the arc from Q to 
Data from the Beaufort Sea, roughly a 500 tan radius about 74°N, lh5°W. 

(From Thorndike and Colony, 19 80. ) 

Strain ellipses from 1979 (solid line) and I960 (broken line) drifting 
buoy data, showing the year long deformation of an initial circle 
(drawn over Greenland). 


(a) The standard deviations of the large scale average velocity deriva- 
tives, as a function of radius of the averaging region, stt ? ) . 

(b) The ratio of the variance of the error in estimating average 
velocity derivatives to the variance in the derivatives ’tVie.wtse.lves, 
as a function of the number of measurements. The dependence on the 
radius of the averaging region is indicated by the shaded width of 
the curves. 

These curves are sensitive to the assumed correlation functions, 
as discussed in the text. 

The total opening and closing versus & . 10 3 realizations of the random 

field with Poisson parameter /\i_, = 8. The velocity discontinuities 
were drawn independently from the unit normal distribution. 

Distributions of the total opening for several ranges of 6 . Data taken 
from figure 24. 
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